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Abstract 



We give a definition of mass for conformally compactifiable initial data 
sets. The asymptotic conditions are compatible with existence of gravita- 
tional radiation, and the compactifications are allowed to be polyhomoge- 
neous. We show that the resulting mass is a geometric invariant, and we 
prove positivity thereof in the case of a spherical conformal infinity. When 
R(g) - or, equivalently, tr g K - tends to a negative constant to order one 
at infinity, the definition is expressed purely in terms of three-dimensional 
or two-dimensional objects. 
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1 Introduction 

In 1958 Trautman [43] (see also [42]) has introduced a notion of energy suitable 
for asymptotically Minkowskian radiating gravitational fields, and proved its 
decay properties; this mass has been further studied by Bondi et al. [9] and 
Sachs [37]. Several other definitions of mass have been given in this setting, 
and to put our results in proper perspective it is convenient to start with a 
general overview of the subject. First, there are at least seven methods for 
defining energy-momentum ("mass" for short) in the current context: 

(i) A definition of Trautman [43], based on the Freud integral [22], that 
involves asymptotically Minkowskian coordinates in space-time. The def- 
inition stems from a Hamiltonian analysis in a fixed global coordinate 
system. 

(ii) A definition of Bondi et al. [9], which uses space-time Bondi coordinates. 

(iii) A definition of Abbott-Deser [1], originally introduced in the context of 
space-times with negative cosmological constant, which (as we will see) 
is closely related to the problem at hand. The Abbott-Deser integrand 
turns out to coincide with the linearisation of the Freud integrand, up to 
a total divergence [20] . 

(iv) The space-time "charge integrals", derived in a geometric Hamiltonian 
framework [13,15,19,27]. A conceptually distinct, but closely related, 
variational approach, has been presented in [10]. 

(v) The "initial data charge integrals", presented below, expressed in terms 
of data (g, K) on an initial data manifold. 

(vi) The Hawking mass and its variations such as the Brown- York mass, using 
two-dimensional spheres. 
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(vii) A purely Riemannian definition, that provides a notion of mass for asymp- 
totically hyperbolic Riemannian metrics [14,44]. 

Each of the above typically comes with several distinct variations. 
Those definitions have the following properties: 

(i) The Bondi mass ms requires in principle a space-time on which Bondi 
coordinates can be introduced. However, a null hypersurface extending to 
future null infinity suffices. Neither of this is directly adapted to an anal- 
ysis in terms of usual spacelike initial data sets, is an invariant under 
Bondi-van der Burg-Metzner-Sachs coordinate transformations. Unique- 
ness is not clear, because there could exist Bondi coordinates which are 
not related to each other by a Bondi-van der Burg-Metzner-Sachs coordi- 
nate transformation. 

(ii) Trautman's definition mx requires existence of a certain class of asymp- 
totically Minkowskian coordinates, with tut being invariant under a class 
of coordinate transformations that arise naturally in this context [43] . The 
definition is obtained by evaluating the Freud integral in Trautman's co- 
ordinates. Asymptotically Minkowskian coordinates associated with the 
Bondi coordinates belong to the Trautman class, and the definition is in- 
variant under a natural class of coordinate transformations. Trautman's 
conditions for existence of mass are less stringent, at least in principle 1 , 
than the Bondi ones, tub equals mx of the associated quasi-Minkowskian 
coordinate system, whenever me exists as well. Uniqueness is not clear, 
because there could exist Trautman coordinates which are not related to 
each other by the coordinate transformations considered by Trautman. 

(iii) The Hawking mass, and its variations, are a priori highly sensitive to the 
way that a family of spheres approaches a cut of J? . (This is a major 
problem one faces when trying to generalise the proof of the Penrose 
inequality to a hyperboloidal setting.) It is known that those masses 
converge to the Trautman-Bondi mass when evaluated on Bondi spheres, 
but this result is useless in a Cauchy data context, as general initial data 
sets will not be collections of Bondi spheres. 

(iv) We will show below that the linearisation of the Freud integral coincides 
with the linearisation of the initial data charge integrals. 

(v) We will show below that the linearisation of the Freud integral does not 
coincide in general with the Freud integral. However, we will also show 
that the resulting numbers coincide when decay conditions, referred to as 

x It is difficult to make a clear cut statement here because existence theorems that lead to 
space-times with Trautman coordinates seem to provide Bondi coordinates as well (though 
perhaps in a form that is weaker than required in the original definition of mass, but still 
compatible with an extension of Bondi's definition). 
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strong decay conditions, are imposed. The strong decay conditions turn 
out to be incompatible with existence of gravitational radiation. 

(vi) We will show below that a version of the Brown- York mass, as well as the 
Hawking mass, evaluated on a specific foliation within the initial data set, 
converges to the Trautman-Bondi mass. 

(vii) We will show below that the Freud integrals coincide with the initial data 
charge integrals, for asymptotically CMC initial data sets on which a 
space-equivalent of Bondi coordinates can be constructed. The proof is 
indirect and uses the result, just mentioned, concerning the Brown- York 
mass. 

(viii) It is important to keep in mind that hyperboloidal initial data sets in gen- 
eral relativity arise in two different contexts: as hyperboloidal hypersur- 
faces in asymptotically Minkowskian space-times on which K approaches 
a multiple of g as one recedes to infinity, or for spacelike hypersurfaces 
in space-times with a negative cosmological constant on which \K\ g ap- 
proaches zero as one recedes to infinity (compare [30] ) . This indicates that 
the Abbott-Deser integrals, which arose in the context of space-times with 
non-zero cosmological constant, could be related to the Trautman-Bondi 
mass. It follows from our analysis below that they do, in fact, coincide 
with the initial data charge integrals under the strong decay conditions. 
In space-times with a cosmological constant, the strong decay conditions 
are satisfied on hypersurfaces which are, roughly speaking, orthogonal to 
high order to the conformal boundary, but will not be satisfied on more 
general hypersurfaces. 

(ix) It has been shown in [14] that the strong decay conditions on g are neces- 
sary for a well defined Riemannian definition of mass and of momentum. 
This appears to be paradoxical at first sight, since the strong decay con- 
ditions are incompatible with gravitational radiation; on the other hand 
one expects the Trautman-Bondi mass to be well defined even if there is 
gravitational radiation. We will show below that the initial data charge 
integrals involve delicate cancelations between g and K, leading to a well 
defined notion of mass of initial data sets without the stringent restrictions 
of the Riemannian definition (which does not involve K). 

(x) Recall that for initial data sets which are asymptotically flat in spacelike 
directions, one can define the mass purely in terms of the induced three- 
dimensional metric. An unexpected consequence of our analysis below 
is therefore that any definition of mass of hyperboloidal initial data sets, 
compatible with initial data sets containing gravitational radiations, must 
involve the extrinsic curvature K in a non-trivial way. 

Let us discuss in some more detail the definition of the Trautman-Bondi 
mass. The standard current understanding of this object is the following: 
one introduces a conformal completion at null infinity of (M, g) and assigns 
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a mass totb to sections of the conformal boundary J? using Bondi coordinates 
or Newman-Penrose coefficients. The Trautman-Bondi mass ttitb has been 
shown to be the unique functional, within an appropriate class, which is non- 
increasing with respect to deformations of the section to the future [16]. A 
formulation of that mass in terms of "quasi-spherical" foliations of null cones 
has been recently given in [5] under rather weak differentiability conditions. 

The above approach raises some significant questions. First, there is a 
potential ambiguity arising from the possibility of existence of non-equivalent 
conformal completions of a Lorentzian space-time (see [12] for an explicit exam- 
ple). From a physical point of view, sections of Scri represent the asymptotic 
properties of a radiating system at a given moment of retarded time. Thus, one 
faces the curious possibility that two different masses could be assigned to the 
same state of the system, at the same retarded time, depending upon which of 
the conformally inequivalent completions one chooses. While there exist some 
partial results in the literature concerning the equivalence question [25,38], no 
uniqueness proof suitable for the problem at hand can be found. The first main 
purpose of this paper is to show that the Trautman-Bondi mass of a section S 
of ^+ is a geometric invariant in a sense which is made precise in Section 5.2 
below. 

Next, recall that several arguments aiming to prove positivity of mTB have 
been given (c/., amongst others, [29,33,36,39]). However, the published proofs 
are either incomplete or wrong, or prove positivity of something which might 
be different from the TB mass, or are not detailed enough to be able to form 
an opinion. The second main purpose of this work is to give a complete proof 
of positivity of rriTB, see Theorems 5.4 and 5.7. 

An interesting property of the Trautman-Bondi mass is that it can be given 
a Hamiltonian interpretation [15]. Now, from a Hamiltonian point of view, it 
is natural to assign a Hamiltonian to an initial data set (^,g, K), where 5^ is 
a three-dimensional manifold, without the need of invoking a four-dimensional 
space-time. If there is an associated conformally completed space-time in which 
the completion 5? of 5? meets ^ + in a sufficiently regular, say differentiable, 
section S, then totb(^) can be defined as the Trautman-Bondi mass of the 
section S. (Hypersurfaces satisfying the above are called hyperboloidal.) It 
is, however, desirable to have a definition which does not involve any space- 
time constructions. For example, for non-vacuum initial data sets an existence 
theorem for an associated space-time might be lacking. Further, the initial 
data might not be sufficiently differentiable to obtain an associated space-time. 
Next, there might be loss of differentiability during evolution which will not 
allow one to perform the space-time constructions needed for the space-time 
definition of mass. Finally, a proof of uniqueness of the definition of mass of 
an initial data set could perhaps be easier to achieve than the space-time one. 
Last but not least, most proofs of positivity use three-dimensional hypersurfaces 
anyway. For all those reasons it seems of interest to obtain a definition of mass, 
momentum, etc., in an initial data setting. The third main purpose of this work 
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is to present such a definition, see (3.11) below. 

The final (3 + l)-dimensional formulae for the Hamiltonian charges turn 
out to be rather complicated. We close this paper by deriving a considerably 
simpler expression for the charges in terms of the geometry of "approximate 
Bondi spheres" near Equations (6.1)-(6.2) below. The expression is similar 
in spirit to that of Hawking and of Brown, Lau and York [10]. It applies to 
mass as well as momentum, angular momentum and centre of mass. 

It should be said that our three-dimensional and two-dimensional versions 
of the definitions do not cover all possible hyperboloidal initial hypersurfaces, 
because of a restrictive assumption on the asymptotic behavior of g tJ Kij, see 
(C.l). This condition arises from the need to reduce the calculational complex- 
ity of our problem; without (C.l) we would probably not have been able to do 
all the calculations involved. We do not expect this condition to be essential, 
and we are planning to attempt to remove it using computer algebra in the fu- 
ture; the current calculations seem pretty much to be the limit of what one can 
calculate by hand with a reasonable degree of confidence in the final formulae. 
However, the results obtained are sufficient to prove positivity of ttt-tb {S) for all 
smooth sections of which bound some smooth complete hypersurface 5? , 
because then 5? can be deformed in space-time to a hypersurface which satisfies 
(C.l), while retaining the same conformal boundary S; compare Theorem 5.7 
below. 

Let us expand on our comments above concerning the Riemannian definition 
of mass: consider a CMC initial data set with tr g K = —3 and A = 0, corre- 
sponding to a hyperboloidal hypersurface in an asymptotically Minkowskian 
space-time as constructed in [2] , so that the g- norm of Kij + gij tends to zero as 
one approaches J? + . It then follows from the vacuum constraint equations that 
R(g) approaches —6 as one recedes to infinity, and one can enquire whether 
the metric satisfies the conditions needed for the Riemannian definition of mass 
for such metrics [14]. Now, one of the requirements in [14] is that the back- 
ground derivatives Vg of g be in L 2 (M). A simple calculation (see Section 4 
below) shows that for smoothly compactifiable \g) this will only be the case 
if the extrinsic curvature x °f the conformally rescaled metric vanishes at the 
conformal boundary d5? . It has been shown in [15, Appendix C.3] that the 
u-derivative of x coincides with the Bondi "news function", and therefore the 
Riemannian definition of mass can not be used for families of hypersurfaces in 
space-times with non-zero flux of Trautman-Bondi energy, yielding an unaccept- 
able restriction. Clearly one needs a definition which would allow less stringent 
conditions than the ones in [14, 19], but this seems incompatible with the ex- 
amples in [19] which show sharpness of the conditions assumed. The answer to 
this apparent paradox turns out to be the following: in contradistinction with 
the asymptotically flat case, in the hyperboloidal one the definition of mass 
does involve the extrinsic curvature tensor K in a non-trivial way. The leading 
behavior of the latter combines with the leading behavior of the metric to give 
a well defined, convergent, geometric invariant. It is only when K (in the A < 
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case) or K + g (in the A = case) vanishes to sufficiently high order that one 
recovers the purely Riemannian definition; however, because the leading order 
of K, or K + g, is coupled to that of g via the constraint equations, one obtains 
- in the purely Riemannian case - more stringent conditions on g — b than those 
arising in the general initial data context. 

We will analyse invariance and finiteness properties of the charge integrals 
in any dimension under asymptotic conditions analogous to those in [14, 19]. 
However, the analysis under boundary conditions appropriate for existence of 
gravitational radiation will only be done in space-time dimension four. 

2 Global charges of initial data sets 
2.1 The charge integrals 

Let g and b be two Riemannian metrics on an n-dimensional manifold M, n > 2, 
and let V be any function there. We set 2 

&ij := 9ij bij • (2-1) 

We denote by D the Levi-Civita connection of b and by Rf the scalar curvature 
of a metric /. In [14] the following identity has been proved: 

v/det^ V(R g - R b ) = di fp<y)) + (* + Q) , (2.2) 

where 

W{V) := 2Vdit^ (Vg^ k g^b m + DWg^ k e jk ) , (2.3) 

s := (-FRic(% + DiDjV - & b Vb ij )g ik gi e e k i , (2.4) 
Q : = V(g^ - W + <? l V^)Ric(% + Q' . (2.5) 

Brackets over a symbol denote anti-symmetrisation, with an appropriate nu- 
merical factor (1/2 in the case of two indices). 3 The symbol Af denotes the 
Laplace operator of a metric /. The result is valid in any dimension n > 2. 
Here Q' denotes an expression which is bilinear in ejj and D^e^, linear in V, 
dV and HessT^, with coefficients which are uniformly bounded in a 6-ON frame, 
as long as g is uniformly equivalent to b. The idea behind the calculation lead- 
ing to (2.2) is to collect all terms in R g that contain second derivatives of the 
metric in <9jU*; in what remains one collects in s the terms which are linear in 

2 The reader is warned that the tensor field e here is not a direct Riemannian counter- 
part of the one in [20]; the latter makes appeal to the contravariant and not the covariant 
representation of the metric tensor. On the other hand dj here coincides with that of [14]. 

3 In general relativity a normalising factor 1/167T, arising from physical considerations, is 
usually thrown in into the definition of IT . From a geometric point view this seems purposeful 
when the boundary at infinity is a round two dimensional sphere; however, for other topologies 
and dimensions, this choice of factor does not seem very useful, and for this reason we do not 
include it in CP. 
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eij, while the remaining terms are collected in Q; one should note that the first 
term at the right-hand-side of (2.5) does indeed not contain any terms linear in 
when Taylor expanded at gij = bij. 

We wish to present a generalisation of this formula — Equation (2.11) below 
- which takes into account the physical extrinsic curvature tensor K and its 
background equivalent K; this requires introducing some notation. For any 
scalar field V and vector field Y we define 

A kl = A kl (b,K) := ££ Y b kl -2VK kl , (2.6) 
A = A(g, K) := V - K kl K kl + (tr g K) 2 ^ -2Y k £>, (k\ - 6 l k ti g K^ 



2Q\X» 
\Jdetg 



=p(9,K) =-J k (g,K)/2 

2G ilv n»X v . (2.7) 



The symbol Jz? denotes a Lie derivative. If we were in a space-time context, then 
Q x fl would be the Einstein tensor density, G^ u would be the Einstein tensor, 
while -n} 1 would be the future directed normal to the initial data hypersurface. 
Finally, an associated space-time vector field X would then be defined as 

X = Vn^ + Y k d k = ^d + (Y k - ^N k )d k , (2.8) 

where iV and iV fc are the lapse and shift functions. However, as far as possible 
we will forget about any space-time structures. It should be pointed out that 
our p here can be interpreted as the energy-density of the matter fields when 
the cosmological constant A vanishes; it is, however, shifted by a constant in 
the general case. 



We set 



P kl := g kl tv g K - K kl , ti g K := g kl K kl , 



with a similar definition relating the background quantities K and P; indices 
on K and P are always moved with g while those on K and P are always moved 
with b. 

We shall say that (V, Y) satisfy the (background) vacuum Killing Initial 
Data (KID) equations if 

A i3 (b,K) = = S kl (b,K) , (2.9) 

where 

S kl = S k \b, K) := V ^2P ml P m k - -^-jtr b PP kl + Ric(6) H - R b b kl \ 

- 5e Y P kl + A b Vb kl - D k D l V . (2.10) 



Vacuum initial data with this property lead to space-times with Killing vectors, 
see [34] (compare [7]). We will, however, not assume at this stage that we are 
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dealing with vacuum initial data sets, and we will do the calculations in a 
general case. 

In Appendix A we derive the following counterpart of Equation (2.2): 



di(W = Jfeig~ [v(p(g,K)-p(b,K)) +s' + Q" 

+Y k ^~g (j k (g,K) - J k (b,K)) , 

(2.11) 

where 

¥(Y) := 2^g~ (P l k - P l k )Y k - ±Y l P mn e mn + l -Y k P l k b mn e mn . 

(2.12) 

Further, Q" contains terms which are quadratic in the deviation of g from b 
and its derivatives, and in the deviations of K from K, while s', obtained by 
collecting all terms linear or linearised in e^, except for those involving p and 
J, reads 



s' = (S kl + B kl )e kl + (P kl -P kl )A kl , 



B kl := - 



ikl pmn a imn % fjkl 



(2.13) 



3 Initial data sets with rapid decay 
3.1 The reference metrics 

Consider a manifold M which contains a region M ext C M together with a 
diffeomorphism 

: M cxt [R, oo) x N , (3.1) 

where is a compact boundaryless manifold. Suppose that on [R, oo) x N we 
are given a Riemannian metric &o of a product form 

b :=-^ + r^h, (3.2) 

as well as a symmetric tensor field Kq; conditions on Kq will be imposed later 
on. We assume that h is a Riemannian metric on N with constant scalar 
curvature R^ equal to 

R= U n -l)(n-2)k, ke{0,±l}, ifn>2, 
h |0, k = l , if n = 2; 

here r is a coordinate running along the [R, oo) factor of [R, oo) x N. Here the 
dimension of A is (n — 1); we will later on specialise to the case n + 1 = 4 but 
we allow a general n in this section. There is some freedom in the choice of k 
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when n = 2, associated with the range of the angular variable ip on N = S , 
and we make the choice k = 1 which corresponds to the usual form of the 
two-dimensional hyperbolic space. When (N, h) is the unit round [n — 1)- 
dimensional sphere (S n_1 , g§n-i), then bo is the hyperbolic metric. 

Pulling-back bo using <J> _1 we define on M ext a reference metric b, 

§*b = b . (3.4) 

Equations (3.2)-(3.3) imply that the scalar curvature Rf, of the metric b is 
constant: 

Rb = Rb = n(n — l)k . 

Moreover, the metric b will be Einstein if and only if h is. We emphasise that 
for all our purposes we only need b on M ext , and we continue b in an arbitrary 
way to M \ M ext whenever required. 

Anticipating, the "charge integrals" will be defined as the integrals of U + V 
over "the boundary at infinity", cf. Proposition 3.2 below. The convergence of 
the integrals there requires appropriate boundary conditions, which are defined 
using the following 6-orthonormal frame {/j}i=i, n on M cxt - 

K l h = re t , * = l,...,n-l, ^ U = V r 2 + k d r , (3.5) 
where the e«'s form an orthonormal frame for the metric h. We moreover set 

9ij ■= 9(fi, fj) , K i:j := K(f u fj) , (3.6) 
etc., and throughout this section only tetrad components will be used. 



3.2 The charges 

We start by introducing a class of boundary conditions for which convergence 
and invariance proofs are particularly simple. We emphasise that the asymp- 
totic conditions of Definition 3.1 are too restrictive for general hypersurfaces 
meeting J? in anti-de Sitter space-time, or - perhaps more annoyingly - for 
general radiating asymptotically flat metrics. We will return to that last case 
in Section 5; this requires considerably more work. 

Definition 3.1 (Strong asymptotic decay conditions). We shall say that the 
initial data (g, K) are strongly asymptotically hyperboloidal if: 



'ijj\ 2 



[ f E (W - <%i 2 + \ K a - ^f) + E \f^9i. 

■ E (\S ij + B^\ 2 + |4f) + E \Jk(j3, K) - J k (b,K)\ 



+ 1 

i,j k 



+ \ P (g,K)-p(b,K)\)ro^dfi g <^, (3.7) 
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3C>0 such that C~ x b{X, X) < g(X, X) < Cb(X, X) . (3.8) 

Of course, for vacuum metrics g and b (with or without cosmological con- 
stant) and for background KIDs (V, Y) (which will be mostly of interest to us) 
all the quantities appearing in the second and third lines of (3.7) vanish. 

For hyperboloids in Minkowski space-time, or for static hypersurfaces in 
anti de Sitter space-time, the V's and Y's associated to the translational Killing 
vectors satisfy 

V = 0(r), y/b#{dV,dV)=0(r) , \Y\ b = 0(r) , (3.9) 

where 6* is the metric on T*M associated to b, and this behavior will be 
assumed in what follows. 

Let Af bo ^ o denote 4 the space of background KIDs: 

K,K ■= Y ) I Ajibo, K ) = = S kl (b , K )} , (3.10) 

compare Equations (2.6) and (2.10), where it is understood that V$ and Y$ 
have to be used instead of V and Y there. The geometric character of (2.6) and 
(2.10) shows that if (Vo,Yq) is a background KID for (bo,Ko), then 

(V := V o$-\Y :=$*Y ) 

will be a background KID for (b = (<$>- r )*b ,K = ($- 1 )*i" )- The intro- 
duction of the (Vo, Yo)'s provides a natural identification of KIDs for different 
backgrounds ((^ybo, i^ 1 )* K ) and 1 )*6 , ($J ^Kq) • We have 

Proposition 3.2. Let the reference metric b on M cxt be of the form (3.4), 
suppose that V and Y satisfy (3.9), and assume that <3? is such that Equations 
(3.7)- (3.8) hold. Then for all (V ,Y ) G N^fa the limits 

H^(V , Y ) := lim / (lP(Vb o S" 1 ) + V{^Y )) dS t (3.11) 

R-^oo J r0 $-i =R 

exist, and are finite. 

The integrals (3.11) will be referred to as the Riemannian charge integrals, 
or simply charge integrals. 

Proof. We work in coordinates on M cxt such that is the identity. For any 
R\ , i?2 we have 

/ (IP + V^dSi = [ (IP + V^dSi + f ^(IP + V*) d n x , (3.12) 

Jr=R 2 Jr=R! J [Ri,R2] xN 

4 We denote by K the background extrinsic curvature on the physical initial data manifold 
M, and by Ko its equivalent in the model manifold [R, oo) x N, Kq '■= &*K. 
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and the result follows from (2.11)-(2.13), together with (3.7)-(3.8) and the 
Cauchy-Schwarz inequality, by passing to the limit R2 — > 00. □ 

In order to continue we need some more restrictions on the extrinsic cur- 
vature tensor K. In the physical applications we have in mind in this section 
the tensor field K will be pure trace, which is certainly compatible with the 
following hypothesis: 

|^-^^k = o(r-"/ 2 ). (3.13) 



Under (3.13) and (3.15) one easily finds from (2.12) that 

lim / Y i {Y)dS i = 2 lim / Vtetb \{P l k - P l k )Y k 

(3.14) 

which gives a slightly simpler expression for the contribution of P to H$. 



Under the conditions of Proposition 3.2, the integrals (3.11) define a linear 
map from N bQ ^ Q to R. Now, each map $ used in (3.4) defines in general a 
different background metric b on M ext , so that the maps H§ are potentially 
dependent 5 upon <I>. (It should be clear that, given a fixed h, (3.11) does 
not depend upon the choice of the frame ti in (3.5).) It turns out that this 
dependence can be controlled: 

Theorem 3.3. Under (3.13), consider two maps & a , a = 1,2, satisfying (3.7) 
together with 

E (\9H ~ <M + \ pi i ~ pi j\) + E \M9ij)\ = 

i,3 i,j,k 

fo(r-™/ 2 ) , if n > 2, 
1 0(r -1 ~ e ) , if n = 2, for some e > 0. 



(3.15) 



Then there exists an isometry A of bo, defined perhaps only for r large enough, 
such that 

H$ 2 {V ,Y ) = (V o A-\ A*Y ) . (3.16) 

Remark 3.4. The examples in [20] show that the decay rate (3.15) is sharp 
when Pij = 0, or when Y l = 0, compare [14]- 

Proof. When K = the result is proved, using a space-time formalism, at the 
beginning of Section 4 in [20]. When Y = this is Theorem 2.3 of [14]. It turns 
out that under (3.13) the calculation reduces to the one in that last theorem, 
and that under the current conditions the integrals of U and V are separately 
covariant, which can be seen as follows: On M ex t we have three pairs of fields: 

(g, K), ((^YM^r 1 )*^) , and ((SJ 1 )***,, (^ 2 _1 )*^o) • 



5 Note that the space of KIDs is fixed, as M bo fc Q is tied to (b ,K ) which are fixed once 
and for all. 
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Pulling back everything by $2 to [R, 00) x N we obtain there 

(($ 2 )*<7, ($2)*^) , ((^r 1 $ 2)*bo, (^r 1 $2)*£o) , and (6 , £ ) • 

Now, ($2)*<7 is simply "the metric g as expressed in the coordinate system 
<3?2'\ similarly for (& 2 )*K, and following the usual physicist's convention we 
will instead write 

(g, K) , (61, £1) := (($r 1 o$ 2 )*6o, (<E>r lo$ 2)*^o) , and (6 2 , #2) = (60, K ) , 

which should be understood in the sense just explained. 

As discussed in more detail in [14, Theorem 2.3], there exists an isometry 
A of the background metric 60 , defined perhaps only for r large enough, such 
that ^j -1 o $ 2 is a composition of A with a map which approaches the identity 
as one approaches the conformal boundary, see (3.21)-(3.22) below. It can be 
checked that the calculation of the proof of [14, Theorem 2.3] remains valid, 
and yields 

H* 2 (V ,0) = H< s>1 (V oA-\0) . (3.17) 

(In [14, Theorem 2.3] Equation (2.10) with Y = has been used. However, 
under the hypothesis (3.13) the supplementary terms involving Y in (3.13) 
cancel out in that calculation.) It follows directly from the definition of H§ 
that 

ff*ioA(0, Y ) = (0, A*Y ) . (3.18) 

Since 

H*(V , Y ) = H*(V , 0) + fT*(0, Y ) , 
we need to show that 

tf$ 2 (0,Y ) = H 9l (0,A*Y ) . (3.19) 

In order to establish (3.19) it remains to show that for all Yq we have 

H* loA {0,Y ) = H* 2 (0,Y ) . (3.20) 

Now, Corollary 3.5 of [20] shows that the pull-back of the metrics by $1 o A has 
the same decay properties as that by $1, so that — replacing 3>i by <J>i o A — to 
prove (3.20) it remains to consider two maps ^j -1 = (r\,v^) and ^ 1 = i r 2, v 2 ) 
(where v A denote abstract local coordinates on N) satisfying 

1 — - 

r 2 =r 1 + o(r 1 2 ), (3.21) 
4 =<+o(r7 (1+f) ), (3.22) 

together with similar derivative bounds. In that case one has, in tetrad com- 
ponents, by elementary calculations, 

ptj-PS^P^-Pjj + oir-"), (3.23) 
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leading immediately to (3.20). We point out that it is essential that P l j appears 
in (3.23) with one index up and one index down. For example, the difference 

P tJ - P tJ = o(r- n / 2 ) , 

transforms as 

Pij - Piij = Pij - P2ij - tr bl A£ c 6i + {r- n ) , 

where 

C=(r 2 -r 1 )^ + B^-^) l} 



A 



□ 



4 Problems with the extrinsic curvature of the con- 
formal boundary 

Consider a vacuum space-time (^#, 4 g), with cosmological constant A = 0, 
which possesses a smooth conformal completion (^#, 4 g) with conformal bound- 
ary J? + . Consider a hypersurface 5? such that its completion 5? in is a 
smooth spacelike hypersurface intersecting transversally, with 5? Pi be- 
ing smooth two-dimensional sphere; no other completeness conditions upon 

or upon ^ + are imposed. In a neighborhood of 5? 'DJ^ + one can introduce 
Bondi coordinates [41], in terms of which 4 g takes the form 

* g = -xVe 2 ? du 2 + 2e 2/3 aT 2 du dx + x~ 2 h A B ( dx A - U A du) ( dx B - U B du) . 

(4.1) 

(The usual radial Bondi coordinate r equals l/x.) One has 

h A B = h AB + x X ab + 0{x 2 ) , (4.2) 

where h is the round unit metric on S 2 , and the whole information about 
gravitational radiation is encoded in the tensor field xab- It has been shown 
in [15, Appendix C.3] that the trace-free part of the extrinsic curvature of 5^V\,f 
within 5? is proportional to %■ I n coordinate systems on 5? of the kind used 
in (3.2) this leads to a 1/r decay of the tensor field e of (2.1), so that the 
decay condition (3.7) is not satisfied. In fact, (2.1) "doubly fails" as the K — K 
contribution also falls off too slowly for convergence of the integral. Thus the 
decay conditions of Definition 3.1 are not suitable for the problem at hand. 

Similarly, let be a space-like hypersurface in a vacuum space-time (^#, g) 
with strictly negative cosmological constant A, with a smooth conformal com- 
pletion (M, 4 g) and conformal boundary J^, as considered e.g. in [21, Section 5]. 
Then a generic smooth deformation of 5? at fixed conformal boundary 5? PI J" 
will lead to induced initial data which will not satisfy (3.7). 
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As already mentioned, in some of the calculations we will not consider the 
most general hypersurfaces compatible with the set-ups just described, because 
the calculations required seem to be too formidable to be performed by hand. 
We will instead impose the following restriction: 

d(ti g K) vanishes on the conformal boundary. (4.3) 

In other words, ti g K is constant on 5? Pi J^, with the transverse derivatives of 
tr g K vanishing there as well. 

Equation (4.3) is certainly a restrictive assumption. We note, however, the 
following: 

(i) It holds for all initial data constructed by the conformal method, both if 
A = [2] and if A < [30], for then d(tr g K) is zero throughout 5?. 

(ii) One immediately sees from the equations in [15, Appendix C.3] that in 
the case A = Equation (4.3) can be achieved by deforming ^ in 
while keeping S^T],^ + fixed, whenever an associated space-time exists. It 
follows that for the proof of positivity of m in space-times with a conformal 
completion it suffices to consider hypersurfaces satisfying (4.3). 

5 Convergence, uniqueness, and positivity of the 
Trautman-Bondi mass 

From now on we assume that the space dimension is three, and that 

A = . 

An extension to higher dimension would require studying Bondi expansions in 
n + 1 > 4, which appears to be quite a tedious undertaking. On the other hand 
the adaptation of our results here to the case A < should be straightforward, 
but we have not attempted such a calculation. 

The metric g of a Riemannian manifold (M, g) will be said to be C k com- 
pactifiable if there exists a compact Riemannian manifold with boundary (M ~ 
MUdooMUdM, 5), where d~M = dMUdooM is the metric boundary of (M,g), 
with dM — the metric boundary of (M,g), together with a diffeomorphism 

i/> : intM -> M 

such that 

= x- 2 g , (5.1) 

where x is a defining function for dooM (i.e., x > 0, {x = 0} = dooM, and dx is 
nowhere vanishing on dooM), with g — a metric which is C k up-to-boundary on 
M. The triple (M, g, x) will then be called a C k conformal completion of (M, g). 
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Clearly the definition allows M to have a usual compact boundary. (M, g) will 
be said to have a conformally compactiftable end M ext if M contains an open 
submanifold M ext (of the same dimension that M) such that (M ext , <?|M ext ) is 
conformally compactifiable, with a connected conformal boundary <9ooM ext . 

In the remainder of this work we shall assume for simplicity that the confor- 
mally rescaled metric g is poly homogeneous and C 1 near the conformal bound- 
ary; this means that g is C 1 up-to-boundary and has an asymptotic expansion 
with smooth expansion coefficients to any desired order in terms of powers of 
x and of lnx. (In particular, smoothly compactifiable metrics belong to the 
polyhomogeneous class; the reader unfamiliar with polyhomogeneous expan- 
sions might wish to assume smoothness throughout.) It should be clear that 
the conditions here can be adapted to metrics which are polyhomogeneous plus 
a weighted Holder or Sobolev lower order term decaying sufficiently fast. In 
fact, a very conservative estimate, obtained by inspection of the calculations 
below, shows that relative 0\ n * x (x A ) error terms introduced in the metric be- 
cause of matter fields or because of sub-leading non-polyhomogeneous behavior 
do not affect the validity of the calculations below, provided the derivatives of 
those error term behave under differentiation in the obvious way (an x deriva- 
tive lowers the powers of x by one, other derivatives preserve the powers). (We 
use the symbol / = ln * x (x p ) to denote the fact that there exists N € N and a 
constant C such that |/| < Cx p (l + | lna^) .) 

An initial data set (M, g, K) will be said to be C k (M) x C e (M) conformally 
compactifiable if (M,g) is C k (M) conformally compactifiable and if K is of the 
form 

K lj = x 3 L ij + t -^-g ij , (5.2) 

with the trace- free tensor in C i (M), and with tr g K in C i (M), strictly 
bounded away from zero on M. We note that (3.15) would have required 
\L\g = o(x 1 / 2 ), while we allow \L\g = 0(1). The slower decay rate is necessary 
in general for compatibility with the constraint equations if the trace-free part 
of the extrinsic curvature tensor of the conformal boundary does not vanish 
(equivalently, if the tensor field \AB in (4.2) does not vanish); this follows from 
the calculations in Appendix C. 

5.1 The Trautman-Bondi four- momentum of asymptotically 
hyperboloidal initial data sets - the four-dimensional 
definition 

The definition (3.11) of global charges requires a background metric b, a back- 
ground extrinsic curvature tensor K, and a map For initial data which are 
vacuum near all these objects will now be defined using Bondi coordinates, 
as follows: Let ( r y,g,K) be a hyperboloidal initial data set, by [17,31] the 
associated vacuum space-time (^, 4 g) has a conformal completion with 
perhaps a rather low degree of differentiability. One expects that (^#, 4 g) will 
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indeed be polyhomogeneous, but such a result has not been established so far. 
However, the analysis of [18] shows that one can formally determine all the 
expansion coefficients of a polyhomogeneous space-time metric on y, as well 
as all their time-derivatives on 5?. This is sufficient to carry out all the cal- 
culations here as if the resulting completion were polyhomogeneous. In all our 
calculations from now on we shall therefore assume that (^#, 4 g) has a poly- 
homogeneous conformal completion, this assumption being understood in the 
sense just explained. 

In 4 g) we can always [18] introduce a Bondi coordinate system (u, x, x A ) 
such that y is given by an equation 

u = a(x,x A ), with a(0,x A ) = , a x (0,x A ) > , (5.3) 

where a is polyhomogeneous. There is exactly a six-parameter family of such 
coordinate systems, parameterised by the Lorentz group (the super translation 
freedom is gotten rid of by requiring that a vanishes on 5? Pi We use the 
Bondi coordinates to define the background 4 b: 

% := -du 2 + 2x~ 2 dudx + x~ 2 h AB dx A dx B . (5.4) 

This Lorentzian background metric 4 6 is independent of the choice of Bondi 
coordinates as above. One then defines the Trautman-Bondi four-momentum 
p M of the asymptotically hyperboloidal initial data set we started with as the 
Trautman-Bondi four-momentum of the cut u = of the resulting J^ + ; recall 
that the latter is defined as follows: Let X be a translational Killing vector of 
4 6, it is shown e.g. in [15, Section 6.1, 6.2 and 6.10] that the integrals 

H(y, X, 4 g, 4 6) := lim / W» X (X, 4 g, %)dS vX 

converge. Here W uX (X, 4 g, 4 b) is given by (B.2). Choosing an ON basis X^ for 
the X's one then sets 

Pli {3>) :=H(y,X^g,%) . 

(The resulting numbers coincide with the Trautman-Bondi four-momentum; we 
emphasise that the whole construction depends upon the use of Bondi coordi- 
nates.) 

5.2 Geometric invariance 

The definition just given involves two arbitrary elements: the first is the choice 
of a conformal completion, the second is that of a Bondi coordinate system. 
While the latter is easily taken care of, the first requires attention. Suppose, for 
example, that a prescribed region of a space-time (^#, g) admits two completely 
unrelated conformal completions, as is the case for the Taub-NUT space-time. 
In such a case the resulting p^s might have nothing to do with each other. 
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Alternatively, suppose that there exist two conformal completions which are 
homeomorphic but not diffeomorphic. Because the objects occurring in the 
definition above require derivatives of various tensor fields, one could a priori 
again obtain different answers. In fact, the construction of the approximate 
Bondi coordinates above requires expansions to rather high order of the metric 
at , which is closely related to high differentiability of the metric at , so 
even if we have two diffeomorphic completions such that the diffeomorphism is 
not smooth enough, we might still end up with unrelated values of p^. 

It turns out that none of the above can happen. The key element of the 
proof is the following result, which is essentially Theorem 6.1 of [14]; the proof 
there was given for C°° completions, but an identical argument applies under 
the hypotheses here: 

Theorem 5.1. Let (M,g) be a Riemannian manifold endowed with two C k , 
k>l and polyhomogeneous conformal compactifications (Mi,g l ,xi) and 
(M 2 ,g 2 ,:r 2 ) with compactifying maps ipi and t/> 2 . Then 

tpi 1 o if) 2 : intM 2 -> intMi 

extends by continuity to a C k and polyhomogeneous conformal up-to-boundary 
diffeomorphism from (M 2 ,g 2 ) t° (-^TbSi); particular M\ and M 2 are dif- 
feomorphic as manifolds with boundary. 

We are ready now to prove definitional uniqueness of four-momentum. Some 
remarks are in order: 

(i) It should be clear that the proof below generalises to matter fields near 
,Jp + which admit a well posed conformal Cauchy problem a la Friedrich, 
e.g. to Einstein- Yang-Mills fields [23]. 

(ii) The differentiability conditions below have been chosen to ensure that 
the conformal Cauchy problem of Friedrich [24] is well posed; we have 
taken a very conservative estimate for the differentiability thresholds, and 
for simplicity we have chosen to present the results in terms of classical 
rather than Sobolev differentiability. One expects that C 1 (M) x C (M) 
and polyhomogeneous CMC initial data (g, L) will lead to existence of a 
polyhomogeneous ^ + ; such a theorem would immediately imply a corre- 
sponding equivalent of Theorem 5.2. 

(iii) It is clear that there exists a purely three-dimensional version of the 
proof below, but we have not attempted to find one; the argument given 
seems to minimise the amount of new calculations needed. Such a three- 
dimensional proof would certainly provide a result under much weaker 
asymptotic conditions concerning both the matter fields and the require- 
ments of differentiability at J^ + . 

Theorem 5.2. Let(M,g,K) be aC 7 (M)xC 6 (M) and polyhomogeneous confor- 
mally compactifiable initial data set which is vacuum near the conformal bound- 
ary, and consider two C (M) and polyhomogeneous compactifications thereof as 
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in Theorem 5.1, with associated four-momenta p^, a = 1,2. Then there exists 
a Lorentz matrix A^" such that 

pi = A/p^ . 

Proof. By the results of Friedrich (see [24] and references therein) the maxi- 
mal globally hyperbolic development (^#, 4 g) of (M, g, K) admits C 4 conformal 
completions {Jt a , A g a ), a = 1,2 with conformal factors O a and diffeomorphisms 
^ a : int Jt ' a — ► Jt such that 

The uniqueness-up-to-conformal-diffeomorphism property of the conformal equa- 
tions of Friedrich together with Theorem 5.1 show that extends by 
continuity to a C 4 -up-to-boundary map from a neighborhood of \|/~ 1 (M) C Jt\ 
to ^#2- Let b a be the Minkowski background metrics constructed near the re- 
spective conformal boundaries as in Section 5.1, we have 

(*i of- 1 )*^ = bi , 

so that o 1 ^>2 1 )* defines a Lorentz transformation between the translational 
Killing vector fields of b\ and 62, and the result follows e.g. from [15, Sec- 
tion 6.9]. □ 



5.3 The Trautman-Bondi four-momentum of asymptotically 

CMC hyperboloidal initial data sets - a three-dimensional 
definition 

Consider a conformally compactifiable initial data set (M, g, K) as defined in 
Section 5, see (5.2). We shall say that (M, g, K) is asymptotically CMC if tr g K 
is in C^M) and if 

the differential of ti g K vanishes on dooM ext • (5.5) 

The vacuum scalar constraint equation (p = in (2.7)) shows that, for C 1 (M) x 
C 1 (M) (or for C 1 (M) x C°(M) and poly homogeneous) conformally compacti- 
fiable initial data sets, Equation (5.5) is equivalent to 

the differential of the Ricci scalar R(g) vanishes on dooM ext . (5.6) 

We wish, now, to show that for asymptotically CMC initial data sets one 
can define a mass in terms of limits (3.11). The construction is closely related 
to that presented in Section 5.1, except that everything will be directly read 
off from the initial data: // a space-time as in Section 5.1 exists, then we 
define the Riemannian background metric b on ,y as the metric induced by 
the metric 4 b of Section 5.1 on the hypersurface u = a, and K is defined as 
the extrinsic curvature tensor, with respect to 4 6, of that hypersurface. The 
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map $ needed in (3.11) is defined to be the identity in the Bondi coordinate 
system above, and the metric 60 is defined to coincide with b in the coordinate 
system above. The four translational Killing vectors of 4 b induce on 
four KIDs (V,Y)^, and one can plug those into (3.11) to obtain a definition of 
four-momentum. However, the question of existence and/or of construction of 
the space-time there is completely circumvented by the fact that the asymptotic 
development of the function a, and that ofh, can be read off directly from g and 
K, using the equations of [15, Appendix C.3]. The method is then to read-off 
the restrictions x\y and x A \,y of the space-time Bondi functions x and x A to 5? 
from the initial data, and henceforth the asymptotic expansions of all relevant 
Bondi quantities in the metric, up to error terms 0\ n * x {x 4 ) (order 0(x 4 ) in the 
smoothly compactifiable case) ; equivalently, one needs an approximation of the 
Bondi coordinate x on 5? up to error terms 0\ n * x (x 5 ). The relevant coefficients 
can thus be recursively read from the initial data by solving a finite number 
of recursive equations. The resulting approximate Bondi function x induces 
a foliation of a neighborhood of the conformal boundary which will be called 
the approximate Bondi foliation. The asymptotic expansion of a provides an 
identification of 5^ with a hypersurface u = a in Minkowski space-time with 
the flat metric (5.4). The Riemannian background metric b is defined to be the 
metric induced by 4 b on this surface, and 3> is defined to be the identity in the 
approximate Bondi coordinates. As already indicated, the KIDs are obtained 
on from the translational Killing vector fields of 4 b. The charge integrals 
(3.11) have to be calculated on the approximate Bondi spheres x = e before 
passing to the limit e — > 0. 

The simplest question one can ask is whether the linearisation of the inte- 
grals (3.11) reproduces the linearisation of the Freud integrals under the pro- 
cedure above. We show in Appendix B that this is indeed the case. We also 
show in that appendix that the linearisation of the Freud integrals does not re- 
produce the Trautman-Bondi four-momentum in general, see Equation (B.21). 
On the other hand that linearisation provides the right expression for p^ when 
the extrinsic curvature x °f the- conformal boundary vanishes. Note that under 
the boundary conditions of Section 3 the extrinsic curvature \ vanishes, and 
the values of the charge integrals coincide with the values of their linearised 
counterparts for translations, so that the calculations in Appendix B prove 
the equality of the 3 + 1 charge integrals and the Freud ones for translational 
background KIDs under the conditions of Section 3. 

The main result of this section is the following: 

Theorem 5.3. Consider an asymptotically CMC initial data set which is C 1 
and poly homogeneously (or smoothly) conformally compactifiable. Let $ be de- 
fined as above and let (V, Y) M be the background KIDs associated to space-time 
translations d^. Then the limits (3.11) H$ ((V,Y) M ) taken along approximate 
Bondi spheres {z = e} C y exist and are finite. Further, the numbers 

p^.= H^((V,Y)^) (5.7) 

coincide with the Trautman-Bondi four-momentum of the associated cut in the 
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e x A 



c x 



Lorentzian space-time, whenever such a space-time exists. 

Proof. We will show that IP + Y x coincides with (6.2) below up to a complete 
divergence and up to lower order terms not contributing in the limit, the result 
follows then from (6.3)-(6.4). It is convenient to rewrite the last two terms in 
(2.12) as 

_ ly*pm ne m n + l -Y k P x k b mn e mn , (5.8) 

so that we can use (E.3)-(E.4) with M = \AB = (3 = N A = there. We further 
need the following expansions (all indices are coordinate ones) 

e k i := b hm (g ml -b ml ), 

e\ = 2(3 + x 3 a x M + O ln * x (x 4 ) , 

= ^ W |C + ^ - x 3 N A - j- 2 xHx CD XCD) U + Om^(x 4 ) , 

= \x 2 a xX AC ^ + a> A - 2x\ x (n a + ^{x CD Xcur) + Om* ,(x 4 ) , 

e A B = x X A b + \x 2 x CD Xcdo A b + x 3 £ A b + O ln * x (x 4 ) , 
with || denoting a covariant derivative with respect to h. One then finds 
-\Y x P mn e mn + l -Y k P x k b mn e mn = Y x ■ O ln * x (x 4 ) + Y B ■ O ln * x (x 5 ) , 
-\Y A P mn e mn + l -Y k P A k b mn e mn = Y x • ln * x (x 5 ) + Y B • O w x (x 4 ) . 

This shows that for Y l which are 0(1) in the {x,x ) coordinates, as is the case 
here (see Appendix F), the terms above multiplied by ^detg = 0(x~ 3 ) will 
give zero contribution in the limit, so that in (2.12) only the first two terms will 
survive. Those are clearly equal to the first two terms in (6.2) when a minus 
sign coming from the change of the orientation of the boundary is taken into 
account. 

On the other hand IP does not coincide with the remaining terms in (6.2), 
instead with some work one finds 



2(Xk - Xk)V - IP = V Vdet^ • 2 g AB g xm D A g mB 



l) Xk + 2 > /d^(V fl - 2 b AB )t\ B 



-2^d^d x Vg^ k e jk 
= vVl(- A x AB \\AB + 0{x 2 )) , (5.9) 

where 2 g AB gBC = $ A C an d similarly for 2 b AB . However, integration over S 2 
yields equality in the limit; here one has to use the fact that V is a linear 
combination of I = and 1 spherical harmonics in the relevant order in x (see 
Appendix F), so that the trace-free part of Vims is 0(x). □ 
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5.4 Positivity 

Our next main result is the proof of positivity of the Trautman-Bondi mass: 

Theorem 5.4. Suppose that (M,g) is geodesically complete without boundary. 
Assume that (M, g, K) contains an end which is C 4 x C 3 , or C 1 and polyho- 
mogeneously, compactifiable and asymptotically CMC. If 

y/gijJW <p£L\M) , (5.10) 
then p^ is timelike future directed or vanishes, in the following sense: 

Po>7E^- (5.11) 

Further, equality holds if and only if there exists a V -covariantly constant spinor 
field on M . 

Remark 5.5. We emphasise that no assumptions about the geometry or the 
behavior of the matter fields except geodesic completeness of (M, g) are made 
on M\M cxt . 

Remark 5.6. In vacuum one expects that equality in (5.11) is possible only 
if the future maximal globally hyperbolic development of (M, g, K) is isomet- 
rically diffeomorphic to a subset of the Minkowski space-time; compare [8] for 
the corresponding statement for initial data which are asymptotically flat in 
space-like directions. When K is pure trace one can use a result of Baum [6] to 
conclude that (M, g) is the three-dimensional hyperbolic space, which implies 
the rigidity result. A corresponding result with a general K is still lacking. 

Before passing to the proof of Theorem 5.4, we note the following variation 
thereof, where no restrictions on ti g K are made: 

Theorem 5.7. Suppose that (M,g) is geodesically complete without boundary. 
If (M,g,K) contains an end which is C 4 x C 3 , or C l and polyhomogeneously, 
compactifiable and which is vacuum near the conformal boundary, then the con- 
clusions of Theorem 5.4 hold. 

Proof. It follows from [15, Eq. (C.84)] and from the results of Friedrich that one 
can deform M near J? + in the maximal globally hyperbolic vacuum develop- 
ment of the initial data there so that the hypotheses of Theorem 5.4 hold. The 
Trautman-Bondi four-momentum of the original hypersurface coincides with 
that of the deformed one by [15, Section 6.1]. □ 

Proof of Theorem 5.4: The proof follows the usual argument as proposed 
by Witten. While the method of proof is well known, there are tedious techni- 
calities which need to be taken care of to make sure that the argument applies. 
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Let D be the covariant-derivative operator of the metric g, let V be the 
covariant-derivative operator of the space-time metric 4 g, and let T) be the 
Dirac operator associated with V along M, 

= yviip 

(summation over space indices only). Recall the Sen-Witten [40,45] identity: 
/ ||v^||2 + (^(p + J l 7,7oW-||^||2= I B'i^dSi, (5.12) 

JM\{r>R) JSr 

where the boundary integrand is 

B\i,) = (V> + yDV,V) 9 • (5-13) 

We have the following: 

Lemma 5.8. Let ip be a Killing spinor for the space-time background metric b, 
and let A^d^ be the associated translational Killing vector field = tp^j^ip. 
We have 



lim 



[ < V, V> + tV'V, V > dSi = -J-Pj.A" . (5.14) 
Js(R) ' 4tt 



Proof. For A^dfj, = do the calculations are carried out in detail in Appendix C. 
For general A^ the result follows then by the well known Lorentz-covariance of 
p^ under changes of Bondi frames (see e.g. [15, Section 6.8] for a proof under 
the current asymptotic conditions). □ 

Lemma 5.9. Let be a spinor field on M which vanishes outside of M cxt , and 
coincides with a Killing spinor for the background metric b for R large enough. 
Then 

Vxb G L 2 (M) . 

Proof. In Appendix D we show that 

Dip £ L? . (5.15) 

We then rewrite (5.12) as 

/ || V ^||2 + ^,( / 9 + J l 7 l 7o)V')= / Wfg+ [ FMdSi. 

JM\{r>R} JM\{r>R} Js R 

(5.16) 

By the dominant energy condition (5.10) the function (tp, {p + J l ji)ip) is non- 
negative. Passing with R to infinity the right-hand-side is bounded by (5.15) 
and by Lemma 5.8. The result follows from the monotone convergence theorem. 

□ 

Lemmata 5.8 and 5.9 are the two elements which are needed to establish 
positivity of the right-hand-side of (5.14) whatever A^, see e.g. [3, 14] for a 
detailed exposition in a related setting, or [4,11,35] in the asymptotically flat 
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context. (For instance, Lemma 5.9 justifies the right-hand-side of the implica- 
tion [14, Equation (4.17)], while Lemma 5.8 replaces all the calculations that 
follow [14, Equation (4.18)]. The remaining arguments in [14] require only 



One also has a version of Theorem 5.4 with trapped boundary, using solu- 
tions of the Dirac equation with the boundary conditions of [26] (compare [28]): 

Theorem 5.10. Let (M,g) be a geodesically complete manifold with compact 
boundary dM, and assume that the remaining hypotheses of Theorem 5.4 or of 
Theorem 5. 7 hold. If dM is either outwards-past trapped, or outwards-future 
trapped, then p^ is timelike future directed: 



6 The mass of approximate Bondi foliations 

The main tool in our analysis so far was the foliation of the asymptotic region 
by spheres arising from space-times Bondi coordinates adapted to the initial 
data surface. Such foliations will be called Bondi foliations. The aim of this 
section is to reformulate our definition of mass as an object directly associated 
to this foliation. The definition below is somewhat similar to that of Brown, 
Lau and York [10], but the normalisation (before passing to the limit) used here 
seems to be different from the one used by those authors. 

Let us start by introducing some notation (for ease of reference we collect 
here all notations, including some which have already been introduced elsewhere 
in the paper): Let 4 g liu be a metric of a spacetime which is asymptotically flat 
in null directions. Let gij be the metric induced on a three-dimensional surface 
y with extrinsic curvature Kij. We denote by 3 r*j£ the Christoffel symbols of 
9ij- 

We suppose that we have a function x > the level sets x = const, of 
which provide a foliation of by two-dimensional submanifolds, denoted by 
5f x , each of them homeomorphic to a sphere. The level set u {x = 0}" (which 
does not exist in y) should be thought of as corresponding to ,y + . The metric 
gij induces a metric 2 gAB on each of these spheres, with area element A = 
V det 2 gAB ■ By Uab we denote the extrinsic curvature of the leaves of the x- 
foliation, kAB = 3 F x ab/V9 XX i w ^ n m ean curvature k = 2 g AB kAB- (The reader 
is warned that in this convention the outwards extrinsic curvature of a sphere in 
a flat Riemannian metric is negative.) There are also the corresponding objects 
for the background (Minkowski) metric 4 b fMU , denoted by letters with a circle. 



trivial modifications.) 



□ 




□ 
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Let X be a field of space-time vectors defined along ,5^ . Such a field can 
be decomposed as X = Y + Vn, where Y is a field tangent to S?, n — the 
unit (n 2 = — 1), future-directed vector normal to . Motivated by (5.9), we 
define the following functional depending upon various objects defined on the 
hypersurface and on a vector field X: 



l - lim / F(X)d 2 x , (6.1) 

DVT X-+0 



167T £-+U / a> 

where 

F(X) = 2 v / det^(P z i - P x ^Y l + 2{Xk - Xk)V. (6.2) 

So far the considerations were rather general, from now on we assume that the 
initial data set contains an asymptotically CMC conformally compactifiable end 
and that x provides an approximate Bondi foliation, as defined in Section 5.3. 
(To avoid ambiguities, we emphasise that we do impose the restrictive condition 
(5.5), which in terms of the function a of (C.2) translates into Equation (C.3).) 
We then have A = A up to terms which do not affect the limit x — > (see 
Appendix E.l) so that the terms containing A above can also be written as 
X(k — k) or X(k — k). We will show that in the limit x — > the functional 
(6.1) will give the Trautman-Bondi mass and momentum, as well as angular 
momentum and centre of mass, for appropriately chosen fields X corresponding 
to the relevant generators of Poincare group. 

To study the convergence of the functional when i-tOwe need to cal- 
culate several objects on S fi . We write the spacetime metric in Bondi-Sachs 
coordinates, as in (C.8) and use the standard expansions for the coefficients 
of the metric (see eg. [15, Equations (5.98)-(5.101)]). The covariant derivative 
operator associated with the metric Kab is denoted by Some intermediate 
results needed in those calculations are presented in Appendix E.l, full details 
of the calculations can be found 6 in [32]. Using formulae (E.2)-(E.5) and the 
decomposition of Minkowski spacetime Killing vectors given in Appendix F we 
get (both here and in Appendix F all indices are coordinate ones): 

T{X time ) = Vl[4M - X CD \\CD + 0{x)\ , 

HXtrans) = ~vVl[4M - X CD \\CD + 0{x)\ , 



T(X rot ) = -VI 



AB ( X A\\C KX XCP)\\A 
£ \ x + 16 



HXboost) = -VI 



+ 6A^4 + ^XAcX CD \\D^j v,B + 0{x) 



X C A\\ C V' A , 1 CD 

h t;X xcdv 



X 



+ (j^(x CD Xcd)\\a + QN A + \xACX CD \\D^j v A + 0{x) 



6 The reader is warned that T there is — T here. 
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The underlined terms in integrands corresponding to boosts and rotations di- 
verge when x tends to zero, but they yield zero when integrated over a sphere. 
In the formulae above we set u — a equal to any constant, except for the last 
one where u — a = 0. Moreover, v is a function on the sphere which is a linear 
combination of £ = 1 spherical harmonics and e AB is an antisymmetric tensor 
(more precise definitions are given in Appendix F). 

In particular we obtain: 

Etb = tJ- lim f F{X time )d 2 x , (6.3) 

167T x^O Jy x 

Ptb = -|- lim I T(X trans )d 2 x , (6.4) 

167T x^O J d y R 

where the momentum is computed for a space-translation generator correspond- 
ing to the function v (see Appendix F). It follows that the integrals of T{X) 
are convergent for all four families of fields X. 



6.1 Polyhomogeneous metrics 

In this section we consider polyhomogeneous metrics. More precisely, we will 
consider metrics of the form (4.1) for which the V, (3, U A and Hab have asymp- 
totic expansions of the form 

Ni 
i j=l 

where the coefficients fij are smooth functions. This means that / can be 
approximated up to terms 0(x N ) (for any N) by a finite sum of terms of the 
form fij(x A )x % log 3 ' x, and that this property is preserved under differentiation 
in the obvious way. 

As mentioned in [15, Section 6.10], allowing a polyhomogeneous expansion 
of Yiab of the form 

fiAB = h A B + xxab + x log xD A b + o(x) 

is not compatible with the Hamiltonian approach presented there because the 
integral defining the symplectic structure diverges. (For such metrics it is still 
possible to define the Trautman-Bondi mass and the momentum [16].) It is also 
noticed in [15, Section 6.10], that when logarithmic terms in hAB are allowed 
only at the x 3 level and higher, then all integrals of interest converge. It is 
therefore natural to study metrics for which hAB has the intermediate behavior 

hAB = h A B + x\AB + x 2 Cab (log x) + 0(x 2 ) , (6.5) 

where Cab (log x) is a polynomial of order N in logx with coefficients being 
smooth, symmetric tensor fields on the sphere. Under (6.5), for the T functional 
we find: 

T(X time ) = Vl[AM - X CD \\CD + 0(x\og N+1 x)] , 
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HXtrans) = ~vVl[AM - X ° C \\ CD + 0(xlog N+1 x)} . 

The only difference with the power-series case is that terms with 0(xlog N+1 x) 
asymptotic behavior appear in the error term, while previously we had 0{x). 
Thus J- reproduces the Trautman-Bondi four-momentum for such metrics as 
well. 

In [15] no detailed analysis of the corresponding expressions for the remain- 
ing generators (i.e., X rot and Xb oost ) was carried out, except for a remark in 
Section 6.10, that the asymptotic behavior (6.5) leads to potentially divergent 
terms in the Freud potential, which might lead to a logarithmic divergence for 
boosts and rotations in the relevant Hamiltonians, which then could cease to 
be well defined. Here we check that the values ^(Xrot) and ^(Xboost) remain 
well defined for certain metrics of the form (6.5). Supposing that h A B is of such 
form one finds (see [15, Appendix C.1.2]) 



QuCab = , 



U A = --x 2 x Ac l]c + x 3 W A (logx) + O(x) 



1 

32 



x 2 X CD XCD + 0(xHog^ x), 



2N 



V = --2M + 0(x\og N+1 x) , 

x 



where W A i}ogx) is a polynomial of order TV + 1 in log a; with coefficients being 
smooth vector fields on the sphere. The PF^'s can be calculated by solving 
Einstein equations (which are presented in convenient form in [18]). 

The calculation in case of polyhomogeneous metrics is very similar to the one 
in case of metrics allowing power-series expansion (see Appendix E.2), leading 
to 



F{X rot ) = -\Th 



.AB 



boost i 



-VI 



X° A\\C 



- 3W A - xW Alx v, B + O(l) 



_ 3Wa _ xWaxx v ,a + 0(1) 
x / 



where lx denotes derivation at constant u. As in the previous section, the 
underlined terms tend to infinity for x — > 0, but give zero when integrated over 
the sphere. However some potentially divergent terms remain, which can be 
rewritten as 



1 



e AB (x 3 W A ) lx v. 



1 



B 



(x 3 W A ) lx v 



To obtain convergence one must therefore have: 

lim f ^e AB (x 3 W A ) lx v B Vh~d 2 x = 
^0Jj?l x 2 
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and 

lim I -^(x 3 W A ) lx v A Vld 2 x = 

for any function v which is a linear combination of the £ = 1 spherical harmon- 
ics. 7 It turns out that in the simplest case N = those integrals are identically- 
zero when vacuum Einstein equations are imposed [18]. 



6.2 The Hawking mass of approximate Bondi spheres 

One of the objects of interest associated to the two dimensional surfaces 5? x is 
their Hawking mass, 



m H (y x 



— (i - — f e-e+d 2 fi 

167T \ 16vr Jc? x 

167T V 16vr Jy x \g xx J \ g xx J J 

We wish to show that for asymptotically CMC initial data the Hawking mass 
of approximate Bondi spheres converges to the Trautman Bondi mass. In fact, 
for a = (a^) let us set 

v(a)(9, ip) = ao + ai— , 

where x % has to be expressed in terms of the spherical coordinates in the usual 
way, and 



with d 2 [i = Xd 2 x. Up to the order needed to calculate the limit of the integral, 
on approximate Bondi spheres satisfying (C.3) we have (see Appendices E.l 
and E.2) 

pXX — . ~ 4x x _ \/K 3 , 



and 



(^Z) 2 = — -^-3x 2 (l-2Mx)-x* X CD llCD + O hl * x (x*), 
V g xx J a )X a )X 



k 2 = JL_W +x 2 {1 _2Mx)+x 3 X CD llCD + O lll * x (x 4 ). 



7 See Appendix F. 
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Therefore 




= -AVI + V / K[2x(4M - X CD \\cd) + ln * x {x 2 )] , 



which, together with 
yields 

p H (a,y x ) = -^ f v(a)[Vl(4M-x CD \\CD) + O ln * x (x)}d 2 x. 

Passing to the limit x — > the X CD \\cd term integrates out to zero, leading to 
equality of the Trautman-Bondi four- momentum Ptb(o>=-^) with the limit of 
PH{a,y x ) defined above. 
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A Proof of (2.11) 

From 

2D l (Y k P l k ) = P kl A kl + Y k J k (g, K) + P kl (2VK kl + &ye kl ) (A.l) 

we get 

VR g + 2D l {Y k P l k ) = Vp(g,K)+Y k J k (g,K) + P kl A kl 

-V (p mn P mn - ^-(tr^P) 2 ) + P kl ^Ye kl 

+Ve kl [2P ml P mk - ^P fcI tr b p) + Q2 , 

(A.2) 
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where, here and below, we use the symbol Qi to denote terms which are 
quadratic or higher order in e and P — P. Moreover, we have 



n Y m 



-e kl J? Y P kl + D m {Y m P kl e kl ) - P kl e kl D m Y m , 
\g k \hi + & Y e kl ) + V(g kl - b kl )K kl + Vtv g K 

1 U 5 V o 

7:9 Ml + T ^bP + Li , 

2 n — 1 



(A.3) 
(A.4) 



2D l {Y k P l k ) = P^A kl + Y k J k (b,K)-2V^P mn P mn --^(tv b P) 2 ^ 



+D l (Y k P l k b mn e r 



l )-b mn e mn D l (Y k P l k ) + Q 1 . (A.5) 



(The term L\ in (A.4) is a linear remainder term which, however, will give a 
quadratic contribution in equations such as (A. 7) below.) Using 



-D x {Y k P\) = V[P mn P mn - 



n — 1 



(tnPf 



l -P^A kl - l -Y k J k {b,K) 



= VR b - Vp(b, K) - ^P kl A kl - h k J k (b, K) 



and (A.2)-(A.5) we are led to 



V(R g -R h ) = v[p{g,K)- P {b,k))+Y k [,J k {g,K)-J k (b,k))- 



(A.6) 



djV 1 
\/detg 



+V 



2P ml P k - 



n — 1 

-c kl .tfyP kl + -e kl 

ykl f}kl 



tT b PP 



kl 



pmn p 



n — 1 



(tr 6 P)< 



ukl 



e k i 



tkl ( r>mn 



pmn \ , \sm j 
r -™-mn i 1 ' J r 



imn a p 



kl 



+(P Kl -P Kl )A M + Q 3 , 
where 

Y l (Y) := 2 Vdet g 



(P l k - P l k )Y k - l -Y l P mn e rnn + l -Y k P l k b mn e ri 



(A.7) 



. (A.8) 



Inserting Equation (2.2) into (A.7) one obtains the following counterpart of 
Equation (2.2) 

9 i (U i (y)+V(y)) = ^g~ [v(p(g,K)-p(b,K)) 

+Y k v ffctg~(j k (g,K)-J k (b,k)^+s' + Q"] , 

(A.9) 

where Q" contains terms which are quadratic in the deviation of g from b and its 
derivatives, and in the deviations of K from K, while s' , obtained by collecting 
all terms linear or linearised in e^-, except for those involving p and J, reads 

s > = (S kl + B kl )e kl + (P kl - P kl )A kl , 
1 



B 



kl ._ 



rjd pmn a rmn a p. 

u r s±m.n u ■ n mn r 



nkl 



(A.10) 
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B 3 + 1 charge integrals vs the Freud integrals 

In this appendix we wish to show that under the boundary conditions of Sec- 
tion 3 the numerical value of the Riemannian charge integrals coincides with 
that of the Hamiltonians derived in a space-time setting in [13, 19]: 

if" = p^ X Q al3 -X"L = d a W^ a + ^^(A)X (S , (B.l) 
W^ A = - detg pa \g a ^8^Xf, a , (B.2) 



87rV l — y P o\y up 
2|det 4 6^| / [detffpoj 7[A u]k 



P ~ levrv+det^^lldet^l 5 9 

= 2 0/37 (y[V ]K ). K 

= 2q^ p % - 2S [ ^gT ~ ld"%V, a , (B.3) 

where a semi-colon denotes the covariant derivative of the metric 4 6, square 
brackets denote anti-symmetrisation (with a factor of 1/2 when two indices are 
involved), as before Qp y = ( K Q aa )~ 1 = Ibirg^/ \f\ det g pa \. Further, ( £ a p(A) is 
the Einstein tensor density eventually shifted by a cosmological constant, 

& a p(A) := ^\detg pa \ ^R a p - \g» v + A£$j (B.4) 

(equal, of course, to the energy-momentum tensor density of the matter fields 
in models with matter, and vanishing in vacuum), while 



P 



+ \Q Xa 9^9a P 5 ap - a , (B.5) 

where by g pi/ we denote the matrix inverse to , and 

L := gT [(K, - B^) (T^ - B* w ) - {T%, - B^) (1^ - S^) + r^] . 

(B.6) 

Finally, is the connection of the background metric, and r pv its Ricci tensor 
- zero in our case. For typesetting reasons in the remainder of this appendix 
we write r\ for 4 6, while the symbol b will be reserved for the space-part of 4 6 
and its inverse. While r\ is flat, the reader should not assume that it takes the 
usual diagonal form. Let us denote by f» the linearisation at A b = rj; we find 

16vrW w ~ (e l m' m - e m j) , (B.7) 

16vrW w fc « (e"V° - e° m ' m ) + eV' - eV°] , (B.8) 

leading to 

lforW^X « (V n 6 im - 6 fc "V n ) e mfc;n , (B.9) 
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16vrW°' fc y fc « (Y l b km - Y k b lm ) (e mk '° - e\, m ) . (B.10) 



Here 



and in the linearised expressions of this appendix all space-time indices are 
raised and lowered with 77. Similarly, 



X^ a 



^(b mn e m J k X ;k -bH mk e mn X°. >k + e m b mk b ln X n . !k ) . (B.ll) 



We also have 



X 1 



i* = J fa - Y l Ki k ) , 



X hk = b lm D k Y m - VK lk 
and we have of course assumed that 



(B.12) 
(B.13) 



X = Vn + Y 



is a background Killing vector field, with Y tangent to the hypersurface of 
interest, so that 

b lm D k Y m + b km b x Y m = 2VK lk . (B.14) 



We recall that 



TV 



V = — - 



0/, 



which gives X° = X., X k = Y k — jjN k , with the lapse and shift given by the 
formulae N = J—. , N k = — \j. We will also need the 3 + 1 decomposition 
of the Christoffel symbols B^ of the four-dimensional background metric in 
terms of those, denoted by F m k i(b), associated with the three-dimensional one: 



B r 



N 



kl 



= T m kl (b) + ^-K lk , B u ok = d k \og N-^-K lk , 



N 



N 



B°- 



1.1 — Klk > ' > /.1 

TV 

Linearising the equation for K \j one finds 



->l. 



D k N l - '—o-D k N - NK l k + ^-N m K mk 
N N 



N 



5K kl := ~N {e\,i + e° /;fc - e kl '°) - ±(N) 2 e 00 K kl , (B.15) 

where the relevant indices have been raised with the space-time background 
metric. This leads to the following formula for the linearised ADM momentum: 

6P l k := 6 l k b mn 5K mn - b ml 6K mk + (b li fc k - S l k K^ e i3 . (B.16) 
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A rather lengthy calculation leads then to 
+ 16ttVW 01 ^ + l6TiW 0l k Y k « 



X 



VbVb^b lm (b t e jrn - b rn e 3i ) + Vbe mn (b mn b lk - b ln b mk )D k V 
+VbD k \Ne° rn {Y l b km - Y k b lm )\ + 2Vb(Y l b km - Y k b lm )5K km 

+Vbe mn VlY k K k m b ln - Y k K k l b mn - Y l K mn = 
VbVb^b lm (b iejm - D m e 3l ) + Vbe mn (b mn b lk - b ln b mk )D k V 
+d k \VbNe° m (Y l b km - Y k b lm )] + 2VbY k 5P l k 

+Vbe mn (P l k Y k b mn - Y l P mn ^j . (B.17) 

The integral over a sphere of (B.17) coincides with the linearisation of the 
integral over a sphere of U* + V 1 , as the difference between those expressions is 
a complete divergence. 

Let us finally show that the linearised expression above does not reproduce 
the Trautman-Bondi mass. It is most convenient to work directly in space-time 
Bondi coordinates (u, x,x A ) rather than in coordinates adapted to 5?. We 
consider a metric of the Bondi form (4.1)-(4.2) with the asymptotic behavior 
(C.11)-(C.13); we have 

e 00 = 1 - xV + x- 2 h AB U A U B = 2Mx + O ln * x (x 2 ) , 

eos = x-\e 2 ? - 1) = 0(1) , 

e A = -x~ 2 h AB U B = ]piA B \\B + 0(x) , 

e A B = x~ 2 (h A B ~ h A B) = x~ x xab + 0(1) , 
e33 = e 3 A = . 

(Recall that h denotes the standard round metric on S 2 .) For translational 
4 6-Killing vector fields X M we have X^ v = 0, hence 



70/ 



Now, taking into account (B.7) and (B.8), the linearised Freud superpotential 
takes the form: 

16^W 03 = ^[(eV-eV)I° + (e°/-eV)I 3 

+ (e 3 A'°-e° A > 3 )X A ] , (B.18) 

The following formulae for the non-vanishing B a p^s for the flat metric (C.5) 
are useful when working out (B.18): 



4 B U AB = X l h A B, A B X AB = xflAB-, 



A B A 



xB 



x o B , 



(B.19) 
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4 T}X 



B xx — 2x 



4 ni 



B A BC = T A Bc(h) . 



(B.20) 



The Killing field corresponding to translations in spacetime can be characterised 
by a function k which is a linear combination of the £ = and 1=1 spherical 
harmonics: 

1 



--x 2 Ak , X A 
2 



xh AB K, B ■ 



With some work one finds the following formula for the linearised superpotential 
for time translations (k = 1) 

16^W 03 = V=v{e A A ' 3 -e 3 A ' A ) 

1 

2' 



= sin# 



4M - - X AB \\ab + X CD duXCD + O ln * x (x) 



The resulting integral reproduces the Trautman-Bondi mass if and only if the 
X CD d u XCD term above gives a zero contribution after being integrated upon; 
in general this will not be the case. 



C Proof of Lemma 5.8 

The object of this appendix is to calculate, for large R, the boundary integrand 
that appears in the integral identity (5.12), for a class of hyperboloidal initial 
data sets made precise in Theorem 5.4. We consider a conformally compactifi- 
able polyhomogeneous initial data set (5f,g,K), such that 

ti g K is constant to second order at dS? . (C.I) 

In (^#, A g) we can always [18] introduce a Bondi coordinate system (u, x, x A ) 
such that y is given by an equation 

u = a(x,x A ), with a(0,x A ) = , a x (0,x A ) > , (C.2) 

where a is polyhomogeneous. It follows from [15, Equation (C.83)] that (C.l) 
is equivalent to 

a,xx\x=o = • (C.3) 

(Throughout this section we will make heavy use of the formulae of [15, Ap- 
pendix C] without necessarily indicating this fact.) Polyhomogeneity implies 
then 

a A = O ln * x (x 3 ) , (C.4) 

where we use the symbol / = 0\ n * x {x p ) to denote the fact that there exists 
iVeN and a constant C such that 

|/| < Cx p {l + \\nx\ N ) . 

We also assume that this behaviour is preserved under differentiation in the 
obvious way: 

\d x f\ < CxP-\l + | lnx\ N ) , \d A f\ < CxP(l + | \nx\ N ) , 
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similarly for higher derivatives. The reason for imposing (C.l) is precisely 
Equation (C.3). For more general a's the expansions below acquire many fur- 
ther terms, and we have not attempted to carry through the (already scary) 
calculations below if (C.3) does not hold. 

Alternatively, one can start from a space-time with a polyhomogeneous J? + 
and choose any space-like hypersurface 5? so that (C.2)-(C.3) hold. Such an 
approach can be used to study the positivity properties of the Trautman-Bondi 
mass, viewed as a function on the set of space-times rather than a function on 
the set of initial data sets. 

We use the Bondi coordinates to define the background 4 fe: 

4 6 = -dn 2 + 2x" 2 dudx + x~ 2 h AB dx A dx B , (C.5) 
so that the components of the inverse metric read 

4 6 «x = x 2 iycx = x 4 A b AB = ^AB 
The metric b induced on ,y takes thus the form 



b = x 



-2 



hAB + O ln * x (x 8 )^ dx A dx B + 2(1 - x 2 a^)a^dxdx A 



+2a 2 ^1 - ^x 2 a x ^ {dx^ 2 



(C6) 



Let e a be a local orthonormal co-frame for the unit round metric h on S 2 
(outside of the south and north pole one can, e.g., use e 1 = d6, e 2 = sinOdip), 
we set 



e a ■= x - x e a , e 3 := x _1 e 3 , 



where 



I (y. 

g 3 := J a jX (2 - x 2 a tX )dx + - ' x a tA dx A . 

v y/a x (2 - x 2 a jX ) 

Assuming (C.3)-(C.4), it follows that the co-frame e 1 is close to being orthonor- 
mal for b: 

b = x ~2 [glgl + g 2g2 + ~3~3 + q w = e l e l + e 2 e 2 + + ln » x (x 4 ) . 

(C.7) 

Here and elsewhere, an equality / = ln * x (x p ) for a tensor field f means that 
the components of f in the coordinates (x,x A ) are 0\ n * x (x p ). 

Recall that in Bondi-Sachs coordinates (u, x, x ) the space-time metric takes 
the form: 

4 g = - x Ve 2 ^du 2 + 2e 2/3 x- 2 dudx + X - 2 h AB (dx A - U A du) (dx B - U B du) . 

(C.8) 
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This leads to the following form of the metric g induced on 5? 



x z g = 2 
+ 
+ 



{h AB U A U B - x 3 Ve 2(3 )a >x a >c + e 2/3 a :C - h CB U B a J dx c dx 
{h AB U A U B - x 3 Ve 2/3 )a >D a >c - 2h CB U B a, D + h CD \ dx c dx D 
{h AB U A U B - x 3 V^) {a x f + 2^a x ] (dx) 2 (C.9) 

Let 7jj be defined as 

x 2 g = j A Bdx A dx B + 2^ xA dx A dx + ^ xx (dx) 2 , 

so that 

1CD ■= h CD ~ 2a :{D h c)B U B + (h AB U A U B - x 3 Ve 2l3 )a, D a, c ■ 
IxC ■= (h AB U A U B - x 3 Ve 2/3 )a x a c + e 2/3 a c - h CB U B a x , 



Ix 



h AB U A U B e~ 2p - x 3 v) a x + 2 



e 2 V 



(CIO) 



If we assume that (y, g, K) is polyhomogeneous and conformally C 1 x C°- 
compactifiable, it follows that 



h AB = h AB (l + —X CD Xcd) + xxab + x 2 ( AB + x 3 £ AB + ln * x (x 4 ) , 

where Q AB and £ AB are polynomials in lnx with coefficients which smoothly 
depend upon the x A, s. By definition of the Bondi coordinates we have det h = 
det h, which implies 



Further, 



h AB X AB = h AB ( AB = 



P = ~^X CD Xcdx 2 + Bx 3 + ln . x {x A ) , 
xV = 1 - 2Mx + O ln * x {x 2 ) , 
h AB U B = --XA B \\ B x 2 + u A x 3 + O ln * x (x 4 ) 



(C.11) 
(C.12) 
(C.13) 



where B and u A are again polynomials in In x with smooth coefficients depend- 
ing upon the x A, s, while || denotes covariant differentiation with respect to the 
metric h. This leads to the following approximate formulae 



(h AB U A U B - = O ln * x (x 5 



\flx~x~ = \[2a 
Let 



1 , 1 



-jXA B \\BX 2 - U A X 3 



h C BU B a jD = Oi n * x (x 5 ) , 
+ O ln * x (x 4 ), (C.14) 



1 " 7( a .* + zX CD Xcd)x 2 + ( ^a x M + B)x 



1 



h ah = h(e a ,e b ) , 



+ O lQ * x (x 4 ) . 
(C.15) 
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where e a is a basis of vectors tangent to S 2 dual to e a , and let fi a b be the 
symmetric root of h a b, 

H a chabH b d = $cd , 



or, in matrix notation, 



t \xh\x = id 



(C.16) 



where t y J c a = \i a c stands for the transpose of /x. Let f a , a = 1,2, be the field of 
local orthonormal co-frames for the metric \iab defined by the formula 



~a ,,a fb 

e — H bj ■ 



We set 



so that 



x 2 9 = f 1 f 1 + f 2 P + f 3 f 3 + O lll * x (x*). 
There exists a matrix M k \ such that 

e k = M k J l . 



Now, 
with 



e 3 = (1 - x z a)f 6 - x l b ai i\f + ln * x (x 



2„\#3 2t ,,a 7b 



1 



a := ~^X CD XCD + [ -a, x M + B)x, 



32 

b a e a = b A dx A := J^- 



la ■ 



t;Xa B \\b-uax 



dx J 



The matrix M is easily calculated to be 

M 3 3 = 1 - x 2 a + ln * x {x A ) , M\ = -x 2 b a ^ a b + Xn * x {x A ) , 
M a 3 = , M a 6 = fi a b . 



Since 

it follows that 



e k = x ~ l e k , f k = x~ x ~f 



e k = M k l f, f l = M k l e k , 



where fi and e& stand for frames dual to /' and e 1 . If we choose 

e 1 = d9 , e 2 = sin #dc^ , 

then 



e 3 := (2 - x 2 a x )dx + 



1 — x a . 



d 



y®~x (2 - x 2 a tX ) 
1 — x 2 a r d 



61 89 a,d a yX (2 — x 2 a yX ) dx 



(C.17) 
(C.18) 
(C.19) 



(C.20) 



(C.21) 
(C.22) 
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1 



e-2 







1 — x a . 



d 



a 



sin# [dip a )X {2 — x 2 a )X ) dx 

1 d 

e3 



We also have the relation 



y/ap (2 - x 2 a iX ) dx 
ek = xe k , f k = xf k ■ 

e, = (M _1 ) fc j/ fc , 



with 



(M~ 1 ) 3 3 



M 3 c 



(M~ 1 ) a 3 = 0, (M _1 ) 3 b 



(C.23) 
(C.24) 



Consider, now, the integrand in (5.12): 

B 3 = - < ip, V 3 ^ + tVV^ >=-<ip, 7 3 7 a VaV' > , 

where the minus sign arises from the fact that we will use a g-orthonormal 
frame /j in which is minus the outer-directed normal to the boundary. Here 
tp is assumed to be the restriction to ,y of a space-time covariantly constant 
spinor with respect to the background metric 4 6: 



It follows that 



vv> = o 



4 (^(/a> /i) " ^(/a, ej)) 7 J 70 + J (^(/a) " ^(/a)) 

This allows us to rewrite i? 3 as 

£ 3 = 1 (#(/«,/,■) - K(f a ,ij)) < V,7 3 7V7oV' > 
-i (c^(/ a ) - ^(/ a )) < V,7 3 7 a f7^ > 
= \ (*(/«, fj) ~ K(f a , ijj) (g^Y a - g^Y 3 ) 

+ \ (W3a(/a) " W3a(/«)) W , 



(C.25) 



(C.26) 



where (W, Y l ) denotes the KID 8 associated to the spinor field ip 

W :=<tp,tp > , Y k =< ^,7 fc 7oV' > • 

8 To avoid a clash of notation with the Bondi function V we are using the symbol W for 
the normal component of the KID here. 
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We use the convention in which 

{y,y} = -2^\ 

with 7 — anti-hermitian, and 7* — hermitian. Because tp is covariantly con- 
stant, W and Y % satisfy the following equation 

diW = K i:j Yi . 

Let e~i be an orthonormal frame for b; we will shortly see that we have the 
following asymptotic behaviors, 

&3a{fa) ~ ^3a(/a) = 0(x 2 ) , 

K(f a J j )-K(f a ,e j ) = 0(x 2 ) , 

Y k = 0{x- 1 ), W = 0{x- 1 ), 

which determines the order to which various objects above have to be expanded 
when calculating B 3 . In particular, these equations show that some non-obvious 
cancelations have to occur for the integral of B 3 to converge. 

Since ej is 6-orthonormal, it holds that 

2uj Hj = b([e k , ei],ij) + b([e k , - b([e h ij],e k ) . (C.27) 

It follows from (C.7) that we can choose e k so that 

h = (Si + O ln * x (x 6 ))e e . (C.28) 
This choice leads to the following commutators: 

[ei, e 2 ] = -x cot 8e 2 + Oi n * x (x 3 ) ei + O ln * x (x 3 )e 2 + ln * x (x A )e 3 , (C.29) 

[e a ,e 3 ] = - [a jX (2 - x 2 a tX )] 1/2 e a - (-^- + x^- + [n * x (x 4 ))e 3 

+0 ln » a; (x 4 )e 1 + ln ^(x 4 )e 2 , (C.30) 
[e,,^-] = Cij k e k , (C.31) 

c 2 i 1 = -ci2 1 =O hl * x (x 3 ) , (C.32) 
c 2 i 2 = -C12 2 = xcot6 + O hl * x {x 3 ) , (C.33) 
C12 3 = -C21 3 = O ln * x (x 4 ) , (C.34) 

C 3 a 3 = -C a3 3 = ^ + ^ + ln * x (x 4 ) , (C.35) 

c 3 a 6 = -ca3 b = [a x (2-x 2 a x )]~ 1/2 5 b a + O ln * x (x i ) . (C.36) 
It follows from (C.19) that we can choose fj so that 

f j = (5 k + O ln * x (x 4 ))f k . (C.37) 
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Let M be the transition matrix from the frame to the frame fj, 

f k = M E k ee , 

if we define the functions dij k as 

[fiJj]=d ij k h, (C.38) 

then we have 

d ml k = MhM i m c ij n {M- x f n + {M- X f 3 f m {Mh) - (M- 1 )* ■■jfiiM* m ) . (C.39) 
Chasing through the definitions one also finds that 

M i j = M i j + O in * x (x 4 ) . (C.40) 

This leads to 

d 3aa = M 3 3 c 3aa + O ln * x (x i ) = 2[a, x (2-x 2 a x )]- 1/2 M 3 3 + O ln * x (x 4 ) . 
Now, 

2^fc/j = Ckij + Ckji — cijk , (C-41) 
2wMj = d/dj + d k ji - dijk , (C.42) 

^Sab = d 3 (ab) - 2^ab3 > (C.43) 
U 3 aa = d 3aa , (C.44) 

ti 3ab = [a A 2 - ^^ x )} _1/2 *a& + O ln * x {x 4 ) = O(l) , (C.45) 

^3a3 = ln . a; (x 3 ) . (C.46) 

Using obvious matrix notation, from (C.16) we obtain 

H = id - | x + + ^ 3 u> + Oi n * ^x 4 ) , (C.47) 

where d and iy are polynomials in In x with coefficients which are x^-dependent 
symmetric matrices. The condition that det fi = 1 leads to the relations 

trx = 0, trd=^trx 2 , tiw = ~ti(xd). 

Using (l-i>~ 1 ) a bf 3 {f J ' b a) = /3(det/x) = 0, the asymptotic expansions (C.33)-(C.36) 
together with (C.23)-(C.24), (C.28) and (C.37) one finds the following contri- 
bution to the first term \ (uj 3a {fa) - W3a(/a)) W in B 3 : 

&3a(/a) - ^3a(/a) = M l a UJ 3a i - L0 3aa 

= M 3 a UJ 3a3 + (> b a - 5 b a )w 3 afe + ^3aa ~ ^3aa + Oln* xfa 4 ) 



X 2 



16 



XABX AB + x ^Ma )X + 2B- itr ( x d)) 



+ ln . a; (x 4 ). (C.48) 
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The underlined term would give a diverging contribution to the integral of B 3 
over the conformal boundary if it did not cancel out with an identical term 
from the K contribution, except when \ = 0- 

We choose now the Killing spinor ijj so that 

d 



Wh + Y = 



dt 



in the usual Minkowskian coordinates, where h is the unit normal to 5^ . This 
leads to 

W = -4== [1 + 0{x 2 )] , Y x = — [1 + 0(x 2 )] , Y A = 0{x 2 ) , 



X 



-1 



+ O(x) , Y a = Q(x) . 



In order to calculate the remaining terms in (C.26) we begin with 

K(f a , fa) ~ K(fa, e a ) = x 2 fi c a fi d a K(e c , e d ) - x 2 fi c a K(e c , e ) + O ln * x (x 4 ) , 

(C.49) 

Using the formulae of [15, Appendix C.3] one finds 

fi c af i d a K(e c ,e d ) = h cd K(e c ,e d ) = h AB K AB + O ln * x (x 2 ) , 



T^ab = T u ab ~ a x T x AB + ln * x (x 3 ) , 
T U AB = x^e' 2 ? (h AB - ^h AB , x ) , 



(C.50) 
(C.51) 



1 



AB 



-2/3 



v - (2%17 B) + d u h AB - 2Vx 2 h AB + Vx 3 h AB:X ) , (C.52) 



where denotes the covariant derivative with respect to the metric h. It 
follows that 



T w AB h AB = x~ l e- 2 ? [2 + xa x (U\ A - 2Vx 2 )] + O lTl * x (x 3 ) 
Further 



(C.53) 



N = 



V2a^ 



l + f3+-x 3 a x V + 0(x 4 ) 



Equation (C.50) specialised to Minkowski metric reads 

T^ab = x' l h AB - xa )X h AB + O hl * x (x 3 ) , 

yielding 



K AB = - 



^j2a )X _ 



hAB ^1 - jX 2 a x } + O lu J,!-'; 



(C.54) 



(C.55) 



We therefore have 

x 2 h AB K AB - x 2 l i h a K ah 



2/3 + \x 2 XCDX CD ~ \x 2 a yX (l - xV) 



+ xa x U A u - — tr ( X d) + O ln * x (x A ) , (C.56) 
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which implies 

(K(f a J a )-K(f a ,e a ))Y 3 = 



2a 



,x 



^XABX AB + 



+ x[2B- 3Ma x + ^a xX AB \\AB ~ ^tr ( X d) 



+ O lQ * x (x 3 ) . 

(C.57) 



Here we have again underlined a potentially divergent term. Using the formulae 
of [15, Appendix C] one further finds 

(K{Lh) - k{Lh)) Y a = o ln * x (x 3 ) , 

which will give a vanishing contribution to B 3 in the limit. Collecting this 
together with (C.57) and (C.48) we finally obtain 



b 3 = \x 2 Um-\ x ab ub ) k; Im ,.(,- :{ ). 



(C.58) 



We have thus shown that the integral of B 3 over the conformal boundary is 
proportional to the Trautman-Bondi mass, as desired. 



D Proof of Lemma 5.9 



From 
we have 



Vrj> = 



-\(K{f t ,fj)-k{f t ,ij)) 



i j 
7 7 70 



We start by showing that 

'uijUO-uJijifi)] 7'7V = C(x 2 ). 
In order to do that, note first 

where 

Aj ki := ujjk(fi) - Ujk(fi) = M l iujj k i - ujjki ■ 



We claim that 



s^ k A jki = 0(x 2 ) , A jkk = 0(x 2 ) 



(D.l) 



(D.2) 
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The intermediate calculations needed for this are as follows: 

[ei,e2] = -xcot6e 2 + O hl * x (x 4 ) (D.3) 
(which follows from (C.29)), 

[e a , e 3 ] = - [2a ^ e a + O hl * x (x 4 ) , (D.4) 

[A, A] = -x cot 6 f 2 + 0(x 2 ) , (D.5) 

[A, f a ] = [2a, x }- 1/2 (f a - \x X b J^ + 0(x 2 ) . (D.6) 
In order to calculate A123 + A312 + A231, we note that 

£l2(/3) +& 3 i(/ 2 ) +^23(/l) = "123 + ^2 ""31a + Ml°"23a + 0(x 2 ) 
= ("123 + ^312 + ^23l) - 7^ (X2 a W 3 l a + Xi a ^23a) + 0{x 2 ) . 

Prom W3 a fe = [2a jX ] -1 ^ 2 5 a b + 0\ n * x (x 3 ) the x terms drops out. Equations (D.3)- 
(D.4) show that 

2u 123 + 2^312 + 2(j 2 3i = b([e 2 , ei], e 3 ) + 6([e 3 , e 2 ], ei) + 6([ei, e 3 ], e 2 ) = 0(x 2 ) . 
Similarly it follows from (D.5)-(D.6) that 

2^23 + 2^312 + 2^231 = g([k A], A) + g([k A], A) + g{[k A], A) = 0(x 2 ) , 

yielding finally 

A123 + A 312 + A 23 i = 0(x 2 ) . 

Next, A 3kk = A 3aa is given by (C.48), and is thus 0{x 2 ). We continue with 
A akk = A abb + A a33 . For a = 1 one finds 

^122 + ^133 = g([fi, A], A) + g([fi, A], A) = -xcot e + o{x 2 ) . 

Further, 

"ll(A)+"l2(A)+"l3(A) = UJn3+Va b Ulab + 0(x 2 ) = UJ lkk - ]<CXa^lab + 0(x 2 ) . 

We have W122 = —xcotO + 0\ n * x (x 3 ) = —xcotO + 0(x 2 ), while wi a b = 0(x 2 ) 
as well, so that the %'s can be absorbed in the error terms, leading to 

Aifcfc = "11(A) + "12(A) + "13(A) - ("122 + "133) = o(a; 2 ) • 
For a = 2 we calculate 

W211 + "233 = <7([A, A], A) + g{[k A], A) = o(x 2 ) , 

and it is easy to check now that A 2kk = 0{x 2 ). This establishes (D.2). 
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To estimate the contribution to (D.l) of the terms involving K we will need 
the following expansions 



d 3 a b 



1 



6°a ~ \x X b a 



v 72 "^ L 2 



(D.7) 



d, 



X 2 b a 



3a 



+ O ln * x {x 3 ) =O hl * x (x 2 



dab 3 = x 2 c ab c b c + O hl * x (x 3 ) = 0{x 2 ) 
Now, it follows from [15, Appendix C.3] that 

tr g K-k(f i ,e i ) = O ln * x (x 2 ) . 

Next, we claim that 

K(fj,e k ) - K(f k ,ij) = O ln * x (x 2 ) . 
We have the following asymptotic formulae for the connection coefficients: 
T" xx = 2x~ l a x + 0(x) , F" xA = 0(x) , 

T^ab = x h AB + -XAB + 0(x) , 
and for the extrinsic curvature tensor: 

K xx = -NT" XX = -x~ 2 [ v / 2q~ + 0(x 2 )] = K xx , 
K xA = -NT» xA = 0(1) , K xA = 0{x) , 



(D.8) 
(D.9) 



K AB = -NT" AB 



x 



-2 



h A B + -xxab + 0(x 2 ) 



K AB = - 



h AB + 0{x 2 



(D.10) 



From the above and [15, Appendix C.3] one obtains the following formulae 

1 



K(e k ,ei) = 7==hl + O in * x (x ) 



(D.ll) 



K(e k ,ei) = -- 



1 



Ski + 7; x Xkl 



a/ 2 "^ [ 2 
where we have set X3k = 0. Further 



+ O lQ * x (x 2 ), (D.12) 



M l k = 5 l k --x X l k + O ln * x (x 2 ) 
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K(fj, e k ) - K(f k , ej) = M l jK{e h e k ) - M l k K{e h ij) 

== (s l j ~ \ x X l ^j % + -y== ($k ~ \ x X l k^j Sij + O ln * x (x 2 ) 



This, together with [15, Equation (C.83)] yields 
tr g K-k(f i ,e i ) = J=+ 



= (s l i - \xx l i} 6u+O hl * x {x 2 ) = O ln * x (x 2 ) . 



We also have 

'K{f h f 3 )-k(f u e 3 ) 



-K{h,h)-K{f h ej)^ 

= -tr g K + K(f k ,e k ) 
L 

+ 7 



so that 
Since 
we obtain 



K(f h fj) - k(f h ej)] 7 V = Om* x (x 2 ) ■ 
y/detg = 0(x 3 ) , < tp, tp >= 0{x' v ) , 



E Asymptotic expansions of objects on ^ 

Throughout this appendix coordinate indices are used. 
E.l Smooth case 

Induced metric g. We write the spacetime metric in Bondi-Sachs coordi- 
nates, as in (C.8), and use the standard expansions for the coefficients of the 
metric (see e.g. [15, Equations (5.98)-(5.101)]). Let 5? be given by uj = const., 
where 

to = u — a{x, x ) . 

Two different coordinate systems will be used: (u, x, x A ) and x, x A ) — co- 
ordinates adapted to ,5^ . To avoid ambiguity two different symbols for partial 
derivatives will be used: the comma stands for the derivative with uj = const, 
and I stands for the derivative with u = const. These two derivatives can be 
transformed into each other: 

A, x = A tx + a lx d u A , 
A jA = A tA + a tA d u A . 

For functions not depending on u {e.g., a) the symbols mean the same. Deriva- 
tions in covariant derivatives ^ and ^ are with to = const. 
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Three-dimensional reciprocal metric. We have the following implicit 
formulae for the three dimensional inverse metric g %3 \ 

9 XA 2 AB 

-9 9xB , 

1 _ g xA 

XjXj I ,-y>,-y> ^ X j\. *} 

AB _ 2 AB , 9 9 

where 2 g AB denotes the matrix inverse to (<?ab)- The calculations get very 
complicated in general. To simplify them we will assume a particular form of 
a, i.e., 

a = const. • x + 0(x 3 ) . (E.l) 

This choice is motivated by the form of a for standard hyperboloid t 2 — r 2 = 1 
in Minkowski spacetime. In that case a = j \/l + x 2 — \ = \x + 0(x 3 ). It is 
further equivalent to the asymptotically CMC condition (4.3). With the above 
assumptions we have the following asymptotic expansions: 

2 /LB 



g AB = X 2 h AB + 0{x 7 ) , 



a , 2 (x AC \\c o t\ta„ {x uu xcd) iia „ \ -Mr if 



,AB„. C 



a' + x a :X 
+0(x A ) , 
1 _ 2a x 



2iV /1 x 



16 



a; 



g 11 x 2 



2 

x cy 

1 + 2/3 ^ + Mx 3 a, + 0(x 4 ) 



gxx 



2a 



,x L 



2 

1 _ 2/3 + - Mx 3 a , + 0(x 4 ) 



1 + /? - -x 2 a ^(1 - 2Mx) + 0(x 4 ) 



4 16 2 x 2 a r v 1 



AB JluAB 



x z h™ + 0(x°) . 



Determinant of the induced metric 2 g. The determinant of the 2-dimensional 
metric induced on the Bondi spheres u = const., x = const., equals 

det (Z2 h ) = det {J 1 ) ■ 
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In our case the u = const., x = const, spheres are not exactly the Bondi ones, 
and one finds 

A := TdeT^ = ^ + <3(x 3 ) , A := VdetV) = ^ + 0(x 6 ) . 
We also note the following purely algebraic identities: 



' • Vdet g , A = ■ Vdet b 



(E.2) 



Lapse and shift. The function N (the lapse) and the vector Si (the shift) 
can be calculated as follows: 

1 



N = 



Hence 



N 



1 



x^f2a~ x 
1 

X' 



9ui 



1 + + jx 2 a >x - ^x 3 a x M + 0(x 4 ) 



S x = ^[l + 2(3- x 2 a /x + 2x 3 a /x M + 0(x 4 )] , 



Sa = \x C A\\c - ^N A - -^x{ X CD Xcd)\\a + 0(x 2 ) 



ChristofFel coeficients. Using the induced metric (and the reciprocal met- 
ric) we can compute the coefficients 3 ^ l jk- The results are: 



— + 7 f ± +Pi x + a x [ -x l M --x + d u p + 0(x 6 ) , 
x 2a iX V 2 2 



1 xA 



-X 



X° ' A\\C ~ ^xN A - ^x(x CD XCD)\\A - 7}XXAC\ ~\\D 



,CD 



2a a ^fliA 2a x a o N 
+ + + — ^ + 0(x 2 ) 



x 2 a ;x 



xa . 



r AB = 7, 1" T7 ^ ~xhAB ~ ~xd u XAB ~ — h O(z) , 



2xa a- 4a 



= ^-xa >x x AC \\c + 0(x 2 ) , 



I^B = -U A B + \x A B + \xa, x d u x A B + \x 2 {^H A B -^^x A B) 

+ ix 2 a,, (xBca u x AC - X AC d uX BC + X AD \\db ~ XBD llDA ) + 0(x 3 ) 



3-pA 



r A BC = r(/i) B c + o(^) 
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Extrinsic curvature. The tensor field Kij can be computed as 

where by Au we denote the covariant derivative of a quantity A with respect 
to gij. These covariant derivatives read: 



%Ct,XX r, „n 1 2 



2a 



+ xj3 lx -2(3- -x a. 



,x 



+ -x 3 a tX M + xa x d u f3 + 0(x 4 ) 



S X \A — T 
x z 



PlA a, A 



a,xA 



+ ^X C A\\C ~ N A 



2x 3 a tX 2x 2 a tX ' 4x 



- ^(x CD Xcd)\\a ~ ^XacX CD \\d + 0{x) , 



S 



A\x 



PlA a,A OL,xA 1 C oat 

+ ~rX a\\c - 3N A 



2x 3 a x 2x 2 a 



Ax' 



3 . pry . 1 s~iy\ 

7^{X XCD)\\A- gXACX 



D 



1 



+ 7;a,xd u (x C a\\c) +0(x) , 



1 



~>A\B 



hAB . XAB 1 I 1 r, . 2\ 

2x^r x + ja- x - r hAB ~ r duXAB + ° {x ] 



Derivatives of with respect to u : 

d u 9xx = 4:X~ 2 a :X d u /3 + 0(x) , 

d u 9xA = ^a x d u (x C ' a\\c) + 0(x) , 

d u 9AB = X^duXAB + 0(1) • 
Hence the extrinsic curvature: 



K x 



\/2a,x 



-. xa xx „ „ 1 2 

-1 h xp )r — a x a x 

2a }X 1 4 



+ 2x z a, x M - xa x d u (3 + 0(x 4 )] , 



K 



xA 



1 c a A ol xA 

^X A\\c + — 4xN A 



\x(x CD Xcd)\\a ~ \xxacX CD \\d + 0(x 2 )] , 



Kab = - 



2x 



h A B {3h AB XAB 3 t 1 2 

h - + ^ xd uXAB + 0(x ) 

xa iX xa.. x la. x 4 2 



132 



The Trautman-Bondi mass 



The trace tr„K = K % i can be easily calculated more precisely as 



1 



ti g K = —[2S\-gVd u . 



or after multiplying by y/det g and changing the covariant divergence to the 
ordinary one: 



After some calculations we get 



tv g K 



/2a- 



3 - 3/? - x(3 lx + 



2a 



- ^x 6 a iX x CD \\CD + xa jX (3 iU + 0(x 4 



ADM momenta. The ADM momenta can be expressed in terms of the 
extrinsic curvature: 

P k i=g ki (guti g K-K u ) . 
Substituting the previously calculated ti g K and Kij we get 

1 



P x , = - 



2-2(3- ^x 2 a x (l - 2Mx) - ^x 3 a xX CD \\ C D + 0(x 4 ) 



px px 

r x r x — 



-2(3 + 3x 6 a x M - ^x 3 a xX CD \\ C D + 0{x 



, (E.3) 



pA 
P B 



+0(x 3 ) , 



(2 - 2(3 - x(3 lx + Y^)5 A B + \x X A b - \x 2 a x d uX A B 

,X 



p a x = -\ x 2 ^r xX AC w + o(x 3 ) , 



—x 



2y/2ax~ 
+Q(x 4 ) , 



X A\\C a,xA 



x 2 a tX 



6N A 



3(x CU Xcd)\\a 1 
16 



CD 

:XACX I ID 



P X A ~ P X A 



—x 



2y/2cT~ x 
+0(x 4 ) 



X° A\\C 



Kx CD Xcd)\\a 1 



6N A 



16 



n d 

-jXACX \\D 



(E.4) 
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The second fundamental form /cab- The extrinsic curvature of the leaves 
of the "2+1 foliation" can be computed from the formula 



kAB 



3px 



AB 



Hence 



kAB 



2x 



1 t 1 

n-AB + 7) XAB 



+ \xh A B ~ \xd u XAB —hAB + 0(x 2 ) 

4 2 xa r 



We need a more accurate expansion of the trace k = 2 g AB kAB- The formula 

kAB 



i _ 9 i AB 

gXX AgLOLO gXX 



can be used. It is convenient to calculate 2 g AB V x ab and 2 g AB T ul ab using the 
expressions for the Christoffel symbols given in [15, Appendix C]: 



^ab = ar _1 e- 2 ^(feAB - ^xh ABl x) - a >x T x A B + 0(x a ) , 



1 



T x ab = --e- 2p {2V{h) (A U B) + d u h AB - 2Vx 2 h AB + Vx 3 h ABlx ) , 

where T>(h)A is a covariant derivative with respect to hAB (we use the (u, x, x A ) 
coordinate system and all the differentiations with respect to x, x A are at 
constant u). Hence: 



2AB r w 



T w ab = 2x- Af3x + x 2 a x (U A u - 2Vx 2 ) + 0{x b ) , 



2 9 AB r* AB = -x 2 [U\ A - 2Vx 2 + 0(x 3 ) 
After substitution of relevant asymptotic expansions: 



k = 



1 



>9 A 



\x A X CD \\CD + - 2Mx) + — - — + 0(* 5 ) 



k = y/2^ 

k — k = \/2a y 



\^X CD \\CD + jx 2 (l - 2Mx) + — - A + (x 4 ) 



(E.5) 
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E.2 The polyhomogenous case 

We give only the most important intermediate results which differ from the 
power-series case: 

g xB = -\x 4 x Ac llc + 1 -x 5 W A + 0(x 5 ), 

S A = \x° A\\c-xW A + 0(x) , 
1 ,C ~, 1 2l 



K xA = ^2cT x \-x C A\\c ~ xW A - -x 2 Wai x + 0(x)\ , 
[x C a\\c ~ 3xW A - x 2 W Alx + 0(x)} . 



F Decomposition of Poincare group vectors into tan- 
gential and normal parts 

The generators of Poincare group are given after [15]: 

X ro t = -£ AB a,AV,Bdu, + e AB v tB d A , 
Xtrans = (~v - xa' A v >A + x 2 va tX )d u) - x 2 vd x + xv A d A , 

Xboost = [xv((a + uj)x + l)a tX - (a + uj)v 

- a' A v t A((a + u)x + l)]d u 

— xv[(a + u)x + l]d x + v' A [(a + to)x + l]d A ■ 



The tensor e AB is defined as 



e AB = ^={A,B}, 
Vh 

where {1, 2} = —{2, 1} = 1 and {1, 1} = {2, 2} = 0. By v we denote a function 
on the sphere which is a combination of I = 1 spherical harmonics. If we 
consider embedding of the sphere into M 3 , then 

v(x ) = v — 

r 

and we get a bijection between functions v and vectors (V) 6 M 3 . 

Let us now decompose a vector field X into parts tangent and normal to 

y : 

X = Y + Vn. 
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Y is a vector tangent to y and n is a unit (n 2 = — 1), future-directed normal 
vector. Setting 

r = 2a 2 - h AB a^a,B , 
we have the following decomposition of respective vectors : 

CD 



» me = ^ r (i- xx a > ca > D + o {x 

^=;(l- XxC °° C ° D +0(x 2 )) , 



time 



a ,A _ xx AC a c _ 



Vrot = -e AB a, A v,B 



CD A 

xx a,ca,Da' 



CD, 



+ 0(x 2 )) , 



2r 



Y? ot = -e AB a A v B • - | 1 - 



1 ( \ XX CD Oi,COi,D 



+ 0(x 2 ) , 



V A _ AB 
Y rot — e V ,B 



1 



CD A 

+ e^acVB --U A - x X AC a c - * X a f> D °" + 0(x 2 ) 



Vtrans — 



X^pT 
1 



( A , x X° D a,ca,D . n( 2 \\ 
I — v — xa : v t A + v ■ V U(x ) I , 



Kan. = ~ (-V ~ Xa^V A + v X ^ Da ^,D + Q{ J 
T \ T 



Y 



trans 



= XV 



- X - (-va> A + vx X AC a c + v ^^£^R a ,A _ xa C V ca A + Q ^ 



V b , 



oast 



— {uj + a)v + xva >x — v' C a ; c(wx + ax + 1) 



+ 



C D 

XX a, cap 
2r 



((lo + a)v + v c a t c) + 0(x 



yx 
boost 



-xv H — 
r 

CD 



(w + q)w + xva >x — v' C a t c(wx + ax + 1) 



+ 



((w + a)w + v' c a, c ) + 0{x 2 ) 



yA 

boost 



v A (tox + ax + 1) 



1 

+ - 

r 



(w + a)va' — xva yX a' + v a^oo (lux + ax + 1) 
X cp a,ca, g a^ + ^ + ^ + + ^ 
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